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Using the force-susceptibility formalism of linear quantum measurements, we study the 
dynamics of signal recycled interferometers, such as LIGO-II. We show that, although the 
antisymmetric mode of motion of the four arm-cavity mirrors is originally described by a 
free mass, when the signal-recycling mirror is added to the interferometer, the radiation- 
pressure force not only disturbs the motion of that "free mass" randomly due to quantum 
fluctuations, but also and more fundamentally, makes it respond to forces as though it were 
connected to a spring with a specific optical-mechanical rigidity. This oscillatory response 
gives rise to a much richer dynamics than previously known for SR interferometers, which 
enhances the possibilities for reshaping the noise curves and, if thermal noise can be pushed 
low enough, enables the standard quantum limit to be beaten. We also show the possibility 
of using servo systems to suppress the instability associated with the optical-mechanical 

Cn . interaction without compromising the sensitivity of the interferometer. 
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Next year a network of broadband ground-based laser interferometers, aimed to detect gravitational waves 
(GWs) in the frequency band 10— 10 4 Hz, will begin operations. This network is composed of the Laser Inter- 
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ferometer Gravitational- wave Observatory (LIGO), VIRGO (whose operation will begin in 2004), GEO 600, 
and TAMA 300 |1|] . Given the anticipated noise spectra and the current estimates of gravitational waves 
from various astrophysical sources |2|, it is plausible but not probable that gravitational waves will be de- 
tected with the first generation of interferometers. The original conception of LIGO included an upgrade of 
LIGO to sensitivities at which it is probable to detect a rich variety of gravitational waves [||. The LIGO 
Scientific Collaboration (LSC) H is currently planning this upgrade to begin in 2006. This second stage 
includes: (i) improvement of the seismic isolation system to push the seismic wall downward in frequency 
to 10 Hz, (ii) improvement of the suspension system to lower the noise in the band between ~ 10 Hz and 
^200 Hz, (hi) increase (decrease) of light power (shot noise) circulating in the arm cavities (^lMWatt), 
(iv) improvement in the optics so that they can handle the increased laser power, and (v) introduction of an 
extra mirror, called a signal-recycling (SR) mirror, at the dark-port output. This upgraded configuration of 
LIGO ("advanced interferometer") is sometimes called LIGO-II and its design is sketched in Fig. pi 

The SR mirror (see Fig. HI) sends the signal coming out the dark port back into the arm cavities; in this 
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FIG. 1. Schematic diagram of a signal recycled interferometer such as LIGO-II. The antisymmetric mode of motion 
of the four arm-cavity mirrors (marked by arrows) is monitored by laser interferometry. A signal-recycling mirror is 
used to feed the signal light back into the arm cavities, while a power-recycling mirror is introduced to feed back into 
the arm cavities the unused laser light coming out the bright port. 



sense it recycles the signal, [j The optical system composed of the SR cavity and the arm cavities forms a 
composite resonant cavity, whose eigenfrequencies and quality factors can be controlled by the position and 
reflectivity of the SR mirror. Near its eigenfrequencies (resonances), the device can gain sensitivity. In fact, 
the initial motivation for introducing the SR cavity was based on the idea of using this feature to reshape 
the noise curves, enabling the interferometer to work either in broadband or in narrowband configurations, 
and improving in this way the observation of specific GW astrophysical sources |2| . Historically, the first 
idea for a narrowband configuration, so-called synchronous or resonant recycling, was due to Drever Ji| and 
was subsequently analyzed by Vinet et al. g. It used a different optical topology from Fig. [y. The original 
idea for the optical topology of Fig. [y was due to Meers [g], who proposed its use for dual recycling - a 
scheme which by recycling the signal light increases the storage time of the signal inside the interferometer 



1 The configuration of LIGO-II will also include a power-recycling (PR) mirror between the laser and the beamsplitter 
(see Fig. hi). This mirror recycles back into the arm cavities the unused laser light coming out the bright port and 
increases the light power at the beamsplitter. Besides this effect, the presence of the PR mirror does not affect the 
derivation of the quantum noise at the dark-port output. Therefore, although in our analysis we assume high light 
power, we do not need to take into account the PR mirror in deducing the interferometer's input-output relation. 



and lowers the shot noise. Later, Mizuno et al. |7]-|9| proposed another scheme called Resonant Sideband 
Extraction (RSE), which also uses the optical topology of Fig. [y but adjusts the SR mirror so that the 
storage time of the signal inside the interferometer decreases while the observation bandwidth increases. In 
general, by choosing appropriate detunings [] of the SR cavity, the optical configuration can be in either of 
the two regimes, or in between. These schemes have been experimentally tested by Freise et al. Jlfj] with the 
30 m laser interferometer in Garching (Germany) , and by Mason |ll|] on a table-top experiment at Caltcch 
(USA). 

All the above mentioned theoretical analyses and experiments of SR interferometers PHll| refer to config- 
urations with low laser power, for which the radiation pressure on the arm-cavity mirrors is negligible and 
the noise spectra are dominated by shot noise. However, when the laser power is increased, the shot noise 
decreases while the effect of radiation-pressure fluctuation increases. LIGO-II has been planned to work at 
a laser power for which the two effects are comparable in the observation band 10-200 Hz B. Therefore, to 
correctly describe the quantum optical noise in LIGO-II, the results so far obtained in the literature |ll| 
must be complemented by a thorough investigation of the influence of the radiation-pressure force on the 
mirror motion. 

Until recently the LIGO-II noise curves were computed using a semiclassical approach 0, which, although 
capable of estimating the shot noise, is unable to take into account correctly the effects of radiation-pressure 
fluctuations. Very recently, building on earlier work of Kimble, Levin, Matsko, Thorne and Vyatchanin 
(KLMTV for short) |§, which describes the initial optical configuration of LIGO/TAMA/VIRGO interfer- 
ometers (so-called conventional interferometers) within a full quantum-mechanical approach, we investigated 
the SR optical configuration (Fig. Q) |l^jlj|. Our analysis revealed important new properties of SR inter- 
ferometers, including: (i) the presence of correlations between shot noise and radiation-pressure noise, (ii) 
the possibility of beating the standard quantum limit (SQL) by a modest amount, roughly a factor of two 
over a bandwidth of A/ ~ / [j and (iii) the presence of instabilities in the optical-mechanical system formed 
by the optical fields and the arm-cavity mirrors. We also noticed J14J that the way the SQL is beaten in 
SR interferometer is quite different from standard quantum-nondemolition (QND) techniques [17]] based on 
building up correlations between shot noise and radiation-pressure noise by (i) injecting squeezed vacuum 
into an interferometer's dark port |18| and/or (ii) introducing two kilometer-long filter cavities into the inter- 
ferometer's output port |19J , |l2|] and applying homodyne detection on the filtered light. Indeed, our analyses 
suggest that the improvement in the noise curves comes largely from the resonant features introduced by 
the SR cavity: whereas the amplitude of the classical output signal is amplified near the resonances, the 
output quantum fluctuation is not strongly affected by them. This way of using resonances to beat the SQL 
was first proposed by Braginsky, Khalili and colleagues in their scheme of "optical bar" GW detectors EG] , 



By detuning of the SR cavity we mean the phase gained by the carrier frequency in the SR cavity, see Sec. [II B 

for details. 
3 This performance refers only to the quantum optical noise. The total noise, which includes also all the other 

sources of noise, such as seismic and thermal noise, can beat the SQL only if thermoelastic noise |lfj can also be 

pushed below the SQL. 



where similarly the test mass is effectively an oscillator whose restoring force is provided by in-cavity optical 
fields. For an "optical bar" the free-mass SQL is irrelevant and we can beat the free-mass SQL using classical 
techniques of position monitoring |20| . 

In Ref. [[Uj our analysis was mainly focused on determining the input-output relations for the electro- 
magnetic quadrature fields in a SR interferometer, and evaluating the corresponding noise spectral density. 
The resonant features of the whole device were discussed only briefly. In the present paper we give a de- 
tailed description of the dynamics of the system formed by the optical fields and the mirrors, we discuss the 
origin of the resonances and their possible instabilities, and we analyze the suppression of the instabilities 
by an appropriate control system. In our analysis we have found the Braginsky-Khalili formalism for linear 



quantum measurements |21| very powerful and intuitive, and we use it throughout this paper. 

This paper is divided into two parts: the formalism and its application. In Sec. O we introduce the 
force-susceptibility formalism and discuss some general features of linear quantum-measurement devices. In 



particular, after briefly commenting in Sec. [I A on general quantum- measurement systems, we derive in 



Sec. II B the equations of motion for linear quantum-measurement devices; in Sec. II C we write down a 
set of conditions on the susceptibilities of linear quantum-measurement systems; in Sec. II D we use these 
conditions to construct an effective description of a quantum- measurement process which allows us to identify 
in a straightforward way the shot noise and the radiation-pressure noise. In the subsequent sections we apply 
the formalism developed in Sec. || to SR interferometers. In Sec. Ill we show that SR interferometers can 
be described by the force-susceptibility formalism and we derive their equations of motion, pointing out the 
existence of a "ponderomotive rigidity" . In Sec. |y| we discuss in detail the oscillatory behavior of the system 
induced by the ponderomotive rigidity, its resonances and instabilities. In Sec. [V] we describe the suppression 



of the instability by a feed-back control system which does not compromise the sensitivity. Finally, Sec. VI 
summarizes our main conclusions. As a foundation for our linear analysis of SR interferometers we summarize 
in Appendix |A| some general properties of linear quantum-mechanical systems. 

II. QUANTUM-MEASUREMENT SYSTEMS 

A. General conditions denning a measurement system 

Following Braginsky and Khalili f21f| , we define a measurement process as a transformation from some 
original classical observable which is unknown, e.g., the gravitational- wave amplitude, into another classical 
observable which is known, e.g., the data stored in the computer. Generally, the system which implements 
this process is composed of a probe V , which is directly coupled to the classical observable to be measured 



(for interferometers this is the antisymmetric mode of motion of the four arm-cavity mirrors, see Sec. ill A), 
and the detector T>, which couples to the probe and produces the output observable (for interferometers 
this is the optical system and the photodetector) . A measurement system is drawn schematically in Fig. |2j. 
Because the probe and the detector are quantum mechanical systems, the overall device is called a quantum- 
measurement device. The output observable Z = S + Q contains a classical part S, which depends on the 
classical observable G to be measured, and some quantum noise Q due to the probe, the detector and their 
mutual interaction. 
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FIG. 2. Schematic diagram of a measurement device. G is the classical observable acting on the probe that we 
want to measure, and Z is the detector's observable which describes the output of the measurement system. 



According to the statistical interpretation of Quantum Mechanics p3[ . the output of a quantum- 
measurement process at different times should be simultaneously measurable. One sufficient condition for 
simultaneous measurability is that the Heisenberg operators of the output observable, Z(t), satisfy n 



Z(h),Z(t 2 ) =0 Vti.ta 



(2.1) 



Henceforth, we shall regard Eq. (2.1) as the condition of simultaneous measurability. Although the condition 



(2.1) was originally introduced by Braginsky et al. |17],|2lJ as the definition of quantum- nondemolition (QND) 
observables (see also Refs. [p4|-p6| ) , we introduce and use it for different purposes, as will become clear in the 
following. If the condition (2.1) is satisfied, then any sample of data i Z(ti), Z(t,2), ■ ■ ■ , Z(t n ) > can be stored 
directly as bits of classical data in a classical storage medium, and any noise from subsequent processing of 
the signal can be made arbitrarily small, i.e. all quantum noises are included in the quantum fluctuations 



of Z(t). We want to discuss the simultaneous measurability condition (2.1) more deeply by pointing out the 
following relation, which was also in part discussed by Unruh |24| and Caves, Thorne, Drever, Sandberg and 
Zimmermann in Sec. IV of Ref. P5| , and reviewed subsequently in Ref. |2q |, although from a different point 
of view. 

Simultaneous-Measurability — Zero-Response Relation: For a Quantum Measurement Device 
(QMD), the simultaneous measurability condition for the output Z(t), i.e. [Z (ti) , Z fo)] = Vti, £2, is 
equivalent to requiring that if the device is coupled to an external system via an interaction Hamiltonian of 
the form V(Z,£) where V is an arbitrary function and £ belongs to the external system, then the back action 
on the QMD does not alter the evolution of the output observable Z . 

Proof of necessity. Q Let us suppose that our QMD with output Z evolves under a Hamiltonian H qmd > 
and that [Z(t), Z(t')\ = for all t, t' . Now let us couple it to an arbitrary external system with Hamiltonian 
Hext via a generic interaction term V{Z , £) as specified above, where £ is an observable of the external 
system. The total Hamiltonian is 



H 



H, 



QMD 



Hex. t )+V(Z,£). 



(2.2) 



We refer to this condition as sufficient since for observables that do not satisfy this condition, there may still exist 

a subspace of the Hilbert space of the system in which these observables are simultaneously measurable. 
5 A similar calculation was carried out by Caves et al. in Sec. IV of Ref. |25| . 



If we treat the two terms in the bracket as the zeroth-order Hamiltonian and the interaction Hamiltonian 
V(Z, £) as a perturbation, by applying the results derived in the Appendix [see Eq. (|A9)] we can write the 
Heisenberg operator of the output variable Z as, 

ft 



^pcrt 

i 

h 



(t) = Z(t) + 



dh 



V(Z(ti),£(ti)),Z(t) 



-' rt i-ti 

dt\ I dt 2 



V(Z(t 2 ),£(t 2 )), V(Z(h),£(h)),Z(t) 



with higher order terms of the form [see Eq. (A£)] 



V(Z(t n ),£(t n )), 



V(Z(t 2 ),£(t 2 )), V(Z(h),£(h)),Z(t) 



(2.3) 



(2.4) 



Here Z(t) and £{t) evolve under the Hamiltonians Hqmd and H^xt, respectively. Because they belong to 
two different Hilbert spaces we have [Z(t), £(t')] = for all t, t' . By assumption, we also have \Z(t\), Z(t 2 )\ = 



Vti, t 2 . Using these two facts, we obtain [V (Z (ti) , £ (ti)) , Z (t 2 )] — Vti, t 2 , and then using Eq. (2.3) 
we derive Z pm t(t) = Z(t). This means that the evolution of Z is not affected by the kind of external coupling 
we introduced. 

Proof of sufficiency. Let us suppose the evolution of Z is not affected by any external system of the form 
specified above. Then, in particular, it must be true for the simple interaction Hamiltonian V(Z , £) — —uZ£^ 
where a is some coupling constant which can vary continuously, e.g., in the interval (0, 1], and we choose a 



classical external coupling £. In this particular case Eq. (2.3) becomes 



WtW = Z{t)-a- 



dtt 



Z(t{),Z{t) f(i!) + 0(a 2 ), 



(2.5) 



with higher order terms of the form: a n [Z(t n ), [• • • , [Z(t 2 ), [Z(ti), Z(t)]] ■■■}]. By assumption the LHS of 
Eq. ( |2.5|) does not change when we vary a. The RHS of Eq. (2^) is a power series in a, and using the 
uniqueness of the Taylor expansion, we deduce that all the terms beyond the zeroth order should vanish 
separately. In particular, the first-order term should vanish and we conclude that [Z(t), Z(t')] = for all 
t,t'. □ 

Let us comment on two interesting aspects of the Simultaneous-Measurability - Zero-Response Relation 
given above. 

• This relation links the abstract quantum mechanical idea of simultaneous measurability to the classical 
dynamics of the measurement device, yielding a simple criterion for the quantum-classical transition: 
the observable which corresponds to the classical output variable should have no response to external 
perturbations directly coupled to it. n We shall use this criterion in our analysis of linear systems in 



Sec. IIC. 



6 By directly coupling to Z we mean the interaction Hamiltonian is of the form V(Z,£), since only this form 
guarantees that Z is the only observable of the device that influences the interaction. 



• This relation is also interesting conceptually. In practice, the result of every measurement is read out 
by coupling the measurement device to another system, and the boundary between the "measurement" 
(still part of the QMD) and the "data analysis" (external to the QMD) occurs at a "stage" at which 
no possible direct coupling to the output observable could change the evolution of the output observ- 
able itself. Otherwise at that stage the "external coupling" should still be considered as part of the 
measurement device. 

Before ending this section, let us compare the point of view followed in this section to the one pursued 
in previous QND analyses p4|-p6[, especially Sec. IV of Ref. ^5|. The authors of Refs. |^5|,^6| followed two 
steps in their discussion. First, they searched for a class of observables A(t) of a quantum-mechanical system 
that can be monitored without adding fundamental noise, deducing a condition for A(t) that coincides with 



Eq. (2.1). They called such observables QND observables. Secondly, they found appropriate interaction 
Hamiltonians describing the coupling between A(t) and a measuring apparatus that do not disturb the 
evolution of A(t) during the measurement process. However, in Refs. pq,Eq| there is no clear distinction 
between what we call the detector and the external measurement system; these two systems are referred 
to together as the measuring apparatus. Thus, the observable A(t) does not necessarily coincide with the 
output Z(t) of our probe-detector system, and for this reason we prefer not to call it a QND observable in 
the sense of Refs. [£§-§6). 

As a final remark, we note that whereas in Refs. pq,Eq| the measuring apparatus and the interaction 
Hamiltonian are indispensable parts of a measurement process, in this paper, by distinguishing the detector 
from the external system, we use the latter only as part of a gedanken experiment, by which we clarify the 
relation between simultaneous measurability and the response to external couplings, which will lead to useful 



properties of linear quantum-measurement devices in Sec. II C 



B. Equations of motion of a linear quantum-measurement system: The force-susceptibility formalism 



Starting in this section we shall focus on linear measurement systems. We shall see in Sec. Ill that GW 
interferometers belong to this class of devices. Our analysis has been inspired by the formalism of linear 
quantum- measurement theory introduced by Braginsky and Khalili (Chaps. V, VI and VII of Ref. gjj) and 
is based on the force-susceptibility description of linearly coupled systems under linearly applied classical 
forces (see, e.g., Sec. 6.4 of Ref. ||). 

In a linear measurement process, the device acts linearly and is linearly coupled to the classical observable 
to be measured (see the Appendix for a precise definition of linear systems). We suppose that the device can 
be artificially divided into two linearly coupled, but otherwise independent subsystems: the probe, which 
is subject to the external classical force we want to measure, and the detector, which yields a classical 
output. More specifically, in our Hamiltonian system the probe is coupled to the external classical force G 
by —yG, where y is some linear observable of the probe, while the probe and the detector are coupled by 
a term — xF, where a; is a generalized (linear) displacement of the probe, and F is a linear observable of 
the detector which describes its back-action force on the probe. In general, the observable x to which the 
external force is coupled and the observable y that the detector directly measures might not be the same. 



G 




FIG. 3. Schematic diagram of a linear measurement system. G is the external classical force acting on the probe 
that we want to measure, x is the linear observable of the probe, F is the linear observable of the detector which 
describes the back-action force on the probe, and Z is the linear observable of the detector which describes the output 
of the overall measurement system. 



However, in our idealized model of GW interferometers (Sec. Ill below), x and y are actually the same 
observable, namely the generalized coordinate of the antisymmetric mode of motion of the four arm-cavity 
mirrors (see Fig. [I] and Sec. Ill A), and F is the radiation-pressure force acting on this mode. Henceforth, 
we shall impose y = x. Finally, we denote by Z the linear observable of the detector which describes the 
output of the entire device. A sketchy representation of the measurement device is drawn in Fig. || The 
linear observables x describing the probe V and Z, F describing the detector T> belong to two different 
Hilbert spaces Tip and TLv, respectively, and the Hilbert space of the combined system is Tip <g) TLt>- The 
Hamiltonian is given by 



H 



Hp-xG) +H V 



xF . 



(2.6) 



We shall now derive the equations of motion of the system composed of the linear observables x, Z and F. 
As a first step in our calculation, we regard the Hamiltonians Hp — xG and Hp, as zeroth order Hamiltonians 
for the subsystems V and T>, respectively, and we treat —xF as a linear coupling between V and T>. Working 
in the Heisenberg picture, we obtain the following equations [see Theorem 4 of the Appendix and Eqs. ( |A12| ), 

dt'C zmFm (t,t')x^(t'), (2.7) 



Z (1 \t) = Z^(t) + - 

h 
n 



?(i) 



(<) 



S(G) 






dt'C F (o )F < 0) {t,t')x {l \t'), 
dt'C x{ G) x{ G){t,t')F {1 \t'). 



(2.8) 
(2.9) 



Here CAB(t,t') is a complex number (C-number), called the (time-domain) susceptibility, and is defined by 
Eq. ( Avllj ) of the Appendix, i.e. 



C AB {t,t')= A(t),B(t') 



(2.10) 



[Henceforth, we shall often use the expressions different-time commutator and time-domain susceptibility 



interchangeably] The superscript (1) in Eqs. (p.7|)-(2.9) denotes time evolution under the total Hamiltonian 
H [Eq. fl2.6D ], the superscript (0) on F(t) and Z(t) denotes time evolution under the free Hamiltonian of the 
detector Htj, while the superscript (G) on x(t) refers to the time evolution under the Hamiltonian H-p — xG, 
which describes the probe under the sole influence of G(t). 



As a second step, we want to relate x^ G '(t) to x^°'(t), which evolves under the free probe Hamiltonian 



H-p. Using Theorem 3 in the Appendix and Eqs. ( A10 ), ( All ), we deduce 



x( G \t)=x^(t) + 



dt C x (o) x (o) (t, t ) G(t ) . 



(2.11) 



Noticing from Eq. ( [2.11 ) that x^ differs from x^ by a time dependent C- number, we get C x (g) x (g) (t, t') = 
C x {o) x {o) (t, t'). Using this fact and inserting Eq. (2.11) into Eq. (2.9), we can relate the Heisenberg operators 



evolving under the full Hamiltonian H to those evolving under the free Hamiltonians of the probe and the 
detector H-p and H-d separately: 

r 

ft, 



zW(t) = z(°»(t) 

j?(D( f ) = j?(0)( t ) 



— oo 

t 



dt'C z(0)Fm (t,t')x^{t'), 

dt'C F{0)F ( 0) (t,t')x^(t'), 



x^(t) =x (0) (t) + 



— oo 

t 



dt' 0,(0^(0) (*, t 1 ) [G(t') + F« (*')] • 



(2.12) 
(2.13) 
(2.14) 



A quantity of special interest for us is the displacement induced on a free probe (without any influence of 



the detector) by G(t), namely the second term on the RHS of Eq. ( 2.11 ). For a GW interferometer this 
displacement is Lh(t), where L is the arm-cavity length and h(t) is the differential strain induced by the 
gravitational wave on the free arm-cavity mirrors (the difference in strain between the two arms). In our 
notation we denote this quantity by 

i /"* 

L h{t) = - / dt' C x (o )x (o) (t, t') G{t') , (2.15) 

" J — oo 

and for a GW interferometer G(t) — (m/A)Lh(t), where to/4 is the reduced mass of the antisymmetric 



mode of motion of the four arm-cavity mirrors (see Sees. Ill A and III B). [Note that each mirror has mass 

TO. ] 

Henceforth, we shall assume that both the probe and the detector have time independent Hamiltonians, 
i.e. both Ht> and H-p are time independent. In this case, as shown in the Appendix, the susceptibilities that 
appear in Eqs. ( 2.12) )— (2.14) depend only on t — t' . By transforming them into the Fourier domain, denoting 
by h(£l) the Fourier transform of h(t) and introducing the Fourier-domain susceptibility 



R AB (n) = -l dre inr C AB (0,-T), 



we derive 



Z (1) (O) = Z (0) (CI) + R ZF {ty x (1) (O) , 

fW(jj) = F<-°\n) + R FF (n)x (1 \Q) , 
x {1 \n) = ^ (0) (n) + l h(n) + R xx (n) f (1) (o) . 

Here and below, to simplify the notation we denote Rzf = Rzi°iF( a ) > Rff = RfWfW) i Rx 



R x (o) 3 



(2.16) 



(2.17) 
(2.18) 
(2.19) 

(o). By 



solving Eqs. (2.17)— (2.19) for the full- evolution operators in terms of the free-evolution ones, we finally get: 



?(1) 



(fi) 



1 



F (1) (tt) = 



1 

i-i? M (n)i? FF (n) 



z( 1 )(«) = zW(f7) 



«**■(«) 



i-i^ x (n)i?ir F (n) 



x (0) (o) + z, h(n) + R xx (n) P (0) (n) 
F {0 \n) + r ff (q) (x {0 \n) + Lh(n) 

£ (0) (fi) + L h{fl) + R xx {Vl) F (0) (fi) 



(2.20) 

(2.21) 
(2.22) 



Let us point out that if the kernel relating the full-evolution operators to the free-evolution ones, i.e. 1/(1 — 
Rxx Rff), contains poles both in the lower and in the upper complex plane [with our definition of Fourier 
transform given by Eq. ( |A14| )] , then by applying the standard inverse Fourier transform to Eqs. ( 2.20J )— 



( 2.22 ), we get that x^'{t), F^'it) and Z^'(t) depend on the gravitational-wave field and the free-evolution 
operators x^°\t), F^°\t) and Z(°\t) both in the past and in the future. However, these are not the correct 
solutions for the real motion. This situation is a very common one in physics and engineering (it occurs for 
example in the theory of linear electronic networks p3] and the theory of plasma waves p7J), and the cure for 
it is well known: in order to obtain the (correct) full-evolution operators x^\t), F^(t) and Z^(t) that only 
depend on the past, we have to alter the integration contour in the inverse- Fourier transform, going above 
all the poles in the complex plane. [In the language of plasma physics we have to use the Landau contours.] 
This procedure, which can be justified rigorously using Laplace transforms [p8|, makes x^ l '(t), F^(t) and 
Z^\t) for many systems (including LIGO-II interferometers) infinitely sensitive to driving forces in the 
infinitely distant past. The reason is simple and well known in other contexts: such quantum-measurement 
systems possess instabilities, which can be deduced from the homogeneous solutions of Eqs. (2.20)-( f2.22| ), 
whose eigenfrequencies are given by the equation 1 — R XX (Q) Rff(VI) = 0. The zeros of the equation 
1 — Rxx(ty Rff(Q) — are generically complex and for unstable systems they have positive imaginary 
parts, corresponding to homogeneous solutions that grow exponentially toward the future. 



C. Conditions defining a linear measurement system in terms of susceptibilities 



As we pointed out in Sec. DIAL in order to be identified as the output of the measurement system, the 
observable Z should satisfy [Z(ti),Z(t2)] = 0, Vti, t%, i.e. the condition of simultaneous measur ability. In 
that section, we have also proved the equivalence between this condition and the condition that any external 
coupling to the measurement system through Z should not change the evolution of Z itself. In the following 
we shall take advantage of this equivalence: By imagining that we couple the linear measurement system to 
some external system through Z and by looking at (possible) changes in Z's evolution, we shall obtain a set 
of conditions for the susceptibilities involving Z . 

Let us first restrict ourselves to the simplest possible external coupling, V = —Z£, where £ is a classical 
external force. The total Hamiltonian ( [2.6|) becomes 



H 



[H V 



xG) +Hx 



xF - Z£ 



(H-p-xG^j + (h-d-ZS) 



xF. 



(2.23) 



To derive the equations of motion for the Hamiltonian (|2.23| ) we apply the procedure used in Sec. [IB to 
deduce the equations of motion for the Hamiltonian ( |2.6| ). First, we consider (H-p — xG) and (Hx> — Z £) 



as zeroth order Hamiltonians and relate the operators Zp Crt , -Fp C rt an d ^pert! which evolve under the full 
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Hamiltonian ( 2.23| ), to the operator x^ G \ which evolves under the Hamiltonian (H-p — xG), and the operators 
Z^ and F( £ \ evolving under the Hamiltonian (H-p — Z £), 



?(i) 



v pcrt 



?(i) 



(*) 



"'pert 



(t) 



-'pertW = -^pertW + T 



— oo 
t 



dt' C F{ c) F ^{t,t')x^l t {t') , 



c (1) 

b pcrt 



{t)=x {G \t) + 



h 



dt' C x ( G ) x i,G)(t,t') F t (t') . 



(2.24) 
(2.25) 
(2.26) 



Second, we relate the operators x^ G ' , Z^ ' and F^ ' to the operators a;'- -', Z^ and F^ which evolve under 
H-p and Hp: 

rt 



z ( f2 t (t) = z^(t) + 

£,©(*) =-F (0) (*) + 

x<°\t)=x< \t) + 



— oo 



dt C Z (o) Z (o){t,t )£(t ) , 
dt C F (o) Z (o){t,t )£(t ) , 

t 

dt C x (o) x (o) (t, t ) G(t ) . 



(2.27) 
(2.28) 
(2.29) 



Noticing that ZpJ. t , -Fp ert and x^ differ from Z™, F^ and x^ only by time dependent C-numbers, 
we obtain the following relations: C z{ £) F (£)(t,t') = C Z ( ) F (o)(t,t'), C F (£) F (£)(t,t') — C F ( ) F (o)(t,t') and 
C xi a) x (G)(t,f) = C x( o) x (o)(t,t'). Then, by inserting Eqs. fl2.27D-fl2.29|) into Eqs. ( |2.24| )-( ^2^) , we deduce the 
equations of motion of Z, F and x under the Hamiltonian ( p. 23 ): 



z%{t) = z^\t) + l -J Jt> 

F^ t (t)=F^(t) + U t dt' 



C zWz m (t, £{t') + C zWpm (t, t') x^> Tt (t') 



C F (o) Z (o) (t,t)£(t) + C F (o) F (o) (t, t ) x t (t ) 



(1) 



41(t) = £ (0) (t) + 1 / dt' C xWx <o> (t, t') G(t') + F™ (f) 



(2.30) 
(2.31) 
(2.32) 



From Eqs. (2.30)-(2.32|) we infer that there are two ways the external force £ can influence the evolution of 



pert" 



?(i 



(i) £ can affect Z^ c [. t directly, through the first term in the bracket of Eq. (2.30), unless C 



(t,t') = 
for all t > t' (and thus for all pairs of t and t'); and (ii) £ can influence the evolution of Zyj Tt indirectly 



'Z(°)Z<°) 

(i) 



affecting the evolution of FpJ. t [first term in the bracket of Eq. (2.31)], and through it the evolution of Xp 0rt 



and ZpJ Tt [second terms in the brackets of Eqs. ( 2.32J ), ( [2.30 )1, unless C F (o) Z (o)(t, if) = for all t > t'. 

Now we are ready to deduce the conditions that must be satisfied in order that the evolution of Z not be 
changed by the external coupling £. In principle the two ways £ affects the evolution of Z may cancel each 
other. However, noticing the fact that case (i) does not depend on the probe (only C Z (o) Z (o) matters), but 
case (ii) does (C x (o) x (o) also matters), we see that the cancellation will not always occur if we assume that, 
whatever probe the detector is coupled to, Z always corresponds to the output of the measurement process. 
Thus both conditions must be satisfied: C Z (o) Z (<>) — and C F {o) Z (a) = 0. 

This argument for both conditions can be made more clear by assigning an "effective mass" fj, to the probe 
and consider a continuous family of probes labeled by \i (for interferometers the family of probes are the 
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family of mirrors with different masses). The susceptibility of the coordinate x depends on the effective mass 

as 

C x (o) x (o) oc — , (2.33) 

which simply says that the probe's response to external forces decreases as its effective mass increases. 
Because Z^ and F^ are operators evolving under the free Hamiltonian of the detector, they do not 
depend on /z. Now consider two cases: First, the limiting case of fi — > oo. Then C x (o) x (o) — ► and from 
Eq. ( [2.32 ) we get x J Tt (t) — x^^t). As a consequence, £ affects the evolution of ZpJ Tt only through the 



first term in the bracket of Eq. (2.30) [see case (i) above], unless C Z (c»z(o) (t, t') = for all pairs of t and t' . 



Second, consider the case of finite mass /x, and then conclude that £ will affect the evolution of Z^J rt only 



through the second term in the bracket of Eq. ( 2.30| ) [see case (ii) above], unless C F (o)z(a)(t, t') = for all 
t > if. 

In conclusion we have found that if, whatever the probe is, Z always corresponds to the output of the 
linear measurement device, then the following conditions must be satisfied 

LQM: f^ )z(0) (MO-[f)W,fHf)]=o vm' 

\ C Fmzm (t, f) =[F(°)(t),Z(°)(t')] =0 VOf. 

In the frequency domain these conditions read 

Rzz(ty=0 = R FZ {n). (2.35) 

It is possible to show that LQM [Eqs. ([2.34)] are also sufficient conditions for the simultaneous measurabil- 



ity condition ( |2.l| ) be satisfied independently of the probe's nature; imagine coupling our linear measurement 
system to an external system with an arbitrary Hamiltonian Hext via a generic coupling V(Z,£), £ being 
an external observable, and check whether the evolution of Z is affected by this coupling. The check can be 
achieved by writing the total Hamiltonian as 



H 



H v -xG) + {H V - Z£ + H_ 



EXT 



xF, (2.36) 



and re-doing all the steps followed earlier in this section. It is helpful to notice that the evolutions of Z 



and F under H-p — Z£ + Hext are the same as those under Hx> 7 once the condition LQM, or Eqs. (2.34) 



is satisfied. The result, after a long calculation is that conditions Q2.34) arc sufficient to guarantee that the 



evolution of Z is unaffected by the coupling. The technical details of the proof are left as an exercise for the 
reader. 

D. Effective description of measurement systems 

It is common to normalize the output observable Z to unit signal — e.g., in the case of GW interferometer, 
it is common to set to unity the coefficient in front of the (classical) observable L h we want to measure so 
the normalized output Z has the form: 

d = £f + Lh 1 (2.37) 
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where TV is the so-called signal-referred quantum noise. The observable O can be easily deduced in the 



frequency domain by renormalizing Eq. ( 2.22 ) 
1 — Rxx(ty Rff(Q) 



O(fi) = 



Rzf(£1) 
Z(°)(Q) 



R ZF (n) 



Rxx{ty 



z (1 \n) 

F^(Q)-R FF (Q) 



g(°)(0) 
RzFip.) 



x {0) (n) + Lh(n) 



that is 



0(0) = z(Q) + fix X (n) j"(n) + £ (0) (0) + l /i(o) . 



(2.38) 



(2.39) 



Here we have introduced two linear observables Z and T defined in the Hilbert space TL-d of the detector, 



Z(0) 



£(°)(n) 



f(tt) = F^iQ,) - R FF (Q) 



Rzf{SI) 
In the time domain the output observable 0(t) reads 

/+oo 
dt'Kit-t^Z^it') 
-OO 

where 

"+°° i - R xx (n) R FF (n) 



zW(n) 

RzfW 



K{t) = f 

J — < 



ifii 



^zf(O) 



do 

2^ 



Thus 



0(t) = Z(t) 



dt'C x m xm (t, t') T{t') + x m (t) + L h(t) 



(2.40) 



(2.41) 



(2.42) 



(2.43) 



Using the two properties given by Eqs. (A17) of the Appendix, and applying the conditions LQM [Eqs. (03)], 
we obtain the following commutation relations for the observables Z(t) and T{t) in the Fourier domain 



= = 



r(Q),j*(rt) 



z(n),P(n') =-2irih5(n-n') 



(2.44) 



or in the time domain: f\ 

Z{t),Z{t') = = T(t),T(t') Vt,t', (2.45) 

Z(t),f{t')\=-ih5{t-t') Vt,t'. (2.46) 

It is interesting to notice that, because the observables Zit) and Tit) satisfy the commutation relations 



(2.45), they can be regarded at different times as describing different degrees of freedom. Moreover, because 



7 Note that, if we use the commutator of Z and T to evaluate the susceptibilities, we find naively that Rrz and 
Rzt are proportional to J dr5(r), which is not a well defined quantity. However, introducing an upper cut-off A 
in the frequency domain we can write 8(t) as 8(r) — sinAr/7rr for A — ► +oo, which is symmetric around the origin. 
With this prescription J dri5(r) = 1/2, and the susceptibilities: Rzz — Rtt = Q,Rtz = 1/2, Rzt — —1/2. 
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of Eq. ( [2.46 ), the observables Z(f) and !F(i) can be seen at each instant of time as the canonical momentum 



and coordinate of different effective monitors (probe-detector measuring devices). Therefore, Z(t) and T{t) 
define an infinite set of effective monitors, indexed by t, similar to the successive independent monitors of 
von Neumann's model |2^| for quantum-measurement processes investigated by Caves, Yuen and Ozawa |29(| . 
However, by contrast with von Neumann's model, the monitors defined by Z(t) and J-(t) at different t's 
are not necessarily independent. They may, in fact, have nontrivial statistical correlations, embodied in the 
relations 

(Z(t) Z{t')) ^ const x S(t - t') , (T{t) T{t')) ^ const x 5{t - t') , (Z(t) £(t')) ^ const x 5(t - t') , (2.47) 

where "( )" denotes the expectation value in the quantum state of the system. These correlations can be 
built up automatically by the internal dynamics of the detector - for example they are present in LIGO-type 
GW interferometers @-0. 



Let us now comment on the origin of the various terms appearing in Eq. (2.43) 



The first term Z(t) describes the quantum fluctuations in the monitors' readout variable [see also 



Eq. (2.4C)] which are independent of the probe. In particular, Z does not depend on the effective mass 
H of the probe. Henceforth, we refer to Z as the effective output fluctuation. For an interferometer, 
the quantum noise embodied in Z is the well-known shot noise. 



• 



The second term in Eq. (2.43]) is the effective response of the output at time t to the monitor's back- 



action force at earlier times t' < t. Since C x (o) x (o) oc l//i this part of the output depends on the 
effective mass of the probe. For GW interferometers the back action is caused by radiation-pressure 
fluctuations acting on the four arm-cavity mirrors. In the following we refer to T as the effective 
back-action or radiation-pressure force. The noise embodied in T is called the back-action noise. [In 
the case of GW interferometers, it is also called the radiation-pressure noise, since the back-action is 
just the radiation-pressure force.] 



• 



The third term in Eq. ([2.43) is the free-evolution operator of the probe's coordinate. In principle, this 



is also a noise term. However, in many cases the free-evolution of the probe coordinate is confined to 
a certain uninteresting frequency range, so if we make measurements outside this range, the noise due 



to the free evolution of the probe will not affect the measurement. We shall see in Sec. [II B that this 
will be the case for GW interferometers, as has been pointed out and discussed at length by Braginsky, 
Gorodetsky, Khalili, Matsko, Thome and Vyatchanin (BGKMTV for short) Q. 



The last term in Eq. ( 2.43| ) is the displacement induced on the probe by the classical observable we 
want to measure. 



Within the effective description of the measurement's renormalized output [Eq. ( 2.43J )], it is instructive to 



analyze how the simultaneous measurability condition [0(ti),0(ta)] = Wi,i2, is enforced by the probe- 
detector interaction. To evaluate explicitly the commutation relations of the observable 0, we notice that in 



Eq. ( 2.43 ) the first two terms always commute with the third term, because they belong to the two different 



Hilbert spaces TC-p and Ti-p- The other terms give 
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0(*i),0(*a) 






^lCulDjIo) (*l>*l) ■^'(*l)) ^(*2) + T 



dt 2 C xi o) x (o) (t 2 , t 2 ) T{t 2 



(2.48) 



Hence, the two-time commutator of O(i) is the sum of two terms: the first term depends solely on detector 
observables, while the second term is just the two-time commutator of the free-probe coordinate x^(t). 
Using the commutation relations of Z(t) and !F(t) given by Eqs. ( 2.45 ), ( 2.46| ) it is straightforward to deduce 
that in Eq. ( [2.48 ) the detector commutator exactly cancels the probe commutator. This clean cancellation is 
a very interesting property of probe-detector kinds of quantum-measurement systems and has been recently 
pointed out and discussed at length by BGKMTV in Ref. §0|. 



III. DYNAMICS OF SIGNAL RECYCLED INTERFEROMETERS: EQUATIONS OF MOTION 

In this section we investigate the dynamics of a SR interferometer, showing that it is a probe-detector 
linear quantum-measurement device as defined and investigated in Sec. ||. 



Laser 




Signal-recycling Q\ 
Mirror 




antisym mode 






X e l 



a, 

FIG. 4. On the left panel we draw a SR interferometer, showing the antisymmetric mode of mirror motion (marked 
by arrows), the dark-port and SR optical fields oh, ■■■, ft and the bright-port heldslji,hi, i — 1,2. The conven- 
tional-interferometer optical scheme is contained inside the dashed box. In the right panel we identify the variables, 
x = Jantisym = (x n i — Xra) — (x e i ~ X e 2), Z and F, describing the dynamics of the SR interferometer. 
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A. Identifying the dynamical variables and their interactions 

In gravitational- wave interferometers composed of equal-length arms (the optical configuration adopted by 
LIGO/VIRGO/GEO/TAMA), laser interferometry is used to monitor the displacement of the antisymmetric 
mode of the four arm-cavity mirrors induced by the passage of a gravitational wave (see Fig. 0). 

Recently Kimble, Levin, Matsko, Thorne and Vyatchanin (KLMTV for short) |L2J described a conventional 
(LIGO-I type) interferometer using a full quantum mechanical approach (see the optical scheme inside the 
dashed box in the left panel of Fig. H). KLMTV [Of showed (as has long been known plj) that in this 
kind of interferometer the antisymmetric mode of motion of the four arm-cavity mirrors and the dark-port 
sideband fields (c; and di [] in Fig. ||) are decoupled from other degrees of freedom, i.e. from other modes 
of motion of the four arm-cavity mirrors and from the bright-port sideband fields (§j and hi in Fig. [|) . As 
a consequence, the dynamics relevant to the output signal and the corresponding noise are described only 
by the antisymmetric mode of motion of the four arm-cavity mirrors and the dark-port sideband fields (see 
Appendix B of KLMTV |1J] for details) . This result remains valid for SR interferometers |l^,|lj| : we only 
need to include in the analysis all the optical fields inside the SR cavity, such as ci , di , ei and fi [but not g~i 
or hi] , and those outside the SR cavity, such as at and bi . 

The coordinate of the antisymmetric mode of motion is defined by KLMTV [see Fig. 3 and Eq. (12) of 
Ref. p2[, and the right panel of Fig. H in our paper] as: 



^antisym 



(X n i-X n2 ) ~ (X e i -X e2 ) , (3.1) 



and we identify it with the dynamical variable x introduced in Sec. [IB [see Eq. (2.9)]. The output of the 
detector can be constructed from two independent output observables, the two quadratures 61 and 62 of 
the outgoing electromagnetic field immediately outside the SR mirror (see the left panel of Fig. [|). If a 
homodyne-detection read-out scheme is implemented, then the output is a linear combination of the two 
quadratures, that is 

&£ = sin C 61 + cos C,b 2 , ( = const , (3-2) 



which is a generic quadrature field. f\ We thus identify the dynamical variable Z introduced in Sec. [IB 



[Eq. (|2.7D 1 as: 

Z C ^\. (3.3) 



In particular, when £ = 7r/2 and £ = we have Z\ = b\ and Z 2 = b 



2- 



8 Here at, bi, ci, . . . with i — 1, 2 stand for the two quadrature operators of the electromagnetic field. This formalism 

was developed by Caves and Schumaker [B2J, adopted by KLMTV |12j and the authors iLjlLJ. 

9 Rigorously speaking, the output is the photocurrent, which in the homodyne detection scheme is almost precisely 

proportional to the output quadrature field, but not quite so; see Ref. BOl and the Appendix of Ref. mM for more 
discussion on this point. 
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The radiation-pressure force acting on the arm-cavity mirrors, and coupled to the antisymmetric mode, 
can be directly related to the dark-port quadrature fields. This result was explicitly derived in Appendix B 
of KLMTV |lj]. As a foundation for subsequent calculations, we shall summarize the main steps of their 
derivation: The force acting on each arm-cavity mirror is 2W/c, where W is the power circulating in each 
arm cavity, which is proportional to the square of the amplitude of the electric field propagating toward the 
mirror. In the arm cavities, the electric field can be decomposed into two parts: the carrier and the sideband 
fields. Introducing the carrier amplitude C and the sideband quadrature operators si,2, we have 

" +0 ° dQ 



E(t) = C COS Wot + COS LUot 



dn 

2^ 



e" 20 *?! 



- sino^ot 



•2tt 



e mt s 2 



(3.4) 



where h.c. stands for Hermitian conjugate. (Note that by writing the carrier field as Ccostuot, we have 
adopted the convention used by KLMTV |L2).) Taking the square of E(t), we obtain 



E 2 (t) = [DC component] + [high frequency component (>ljo) 



C 



2tt 



e * m si 



h.c. 



(quadratic terms in si, S2) ■ 



(3.5) 



where we have used the fact that in the integral £1 < u)q. The DC and u)$ ~ 10 15 sec l components are not in 
the detection band of GW interferometers, 10 Hz < J7/27T < 10 4 Hz; in practice they will be counteracted by 



control systems. We also ignore the quadratic terms in Eq. (3.5), since they are much smaller than the linear 
terms. Thus, modulo a factor of proportionality, we obtain in the Fourier domain the following expression 
for the radiation-pressure force acting on each mirror: 



F RP (fi) ex Csi(ft) 



(3.6) 



As shown in Appendix B of Ref. ]12J] , the in-cavity quadrature field si is a combination of the incoming 
quadratures from both the dark and the bright ports. However, the contribution from the bright-port fields 
do not couple to the antisymmetric mode, so the force acting on the antisymmetric mode is due only to 
the incoming fields from the dark port. More specifically, in Sec. 4 of Appendix B of Ref. |12], KLMTV 
related the in-cavity carrier amplitude C and the sideband quadrature si (which they denoted by j\ q|) 
to the input carrier amplitude and ingoing dark-port quadrature C\ (which they denoted by a\). Although 
they did not give the explicit expression we need here for -Frp, it is straightforward to recover it. Using the 
arrows indicated in the right panel of Fig. as positive directions, we findF] 



Frp 



2I hu>o il3 ^ 
■e lp ci, 



(f7 2 + 7 2 )L 2 



(3.7) 



10 We ignore the effect of the arm-cavity optical losses, thus in this case the quadratures ji and hi in Ref. |12| are 

equal. 

11 This result can be obtained from Eq. (B21) of KLMTV jjjj using the fact that x B a = —l/mQ. 2 F RP . Since in 

this paper we ignore optical losses, in Eq. (B21) we can replace /3* and A3* by /3 and /C and ignore the noise operator 
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where Lj{) is the carrier laser frequency, Iq is the carrier light power entering the beamsplitter, 2/3 = 
2 arctan fi/7 is the net phase gained by the sideband frequency tt while in the arm cavity, 7 = Tc/AL 
is the half bandwidth of the arm cavity (T is the power transmissivity of the input mirrors and L is the 



length of the arm cavity). We identify the force F-^p with the dynamical variable F introduced in Sec. [IB 



[see Eq. ( [2.81) 1 



Applying Newton's law to the four mirrors, we deduce 

m x = AF + other forces , (3-9) 

where "other forces" refer to forces not due to the optical-mechanical interaction, e.g., the force due to the 
gravitational wave and thermal forces. By identifying the reduced mass of the antisymmetric mode as m/4, 



we obtain that the coupling term in the total Hamiltonian (2.6) is —xF. [The reduced mass coincides with 
the effective mass of the probe /x introduced in Sec. II.] 

Note that, by assuming the four forces acting on the arm-cavity mirrors are equal, we have made the 
approximation used by KLMTV [O] of disregarding the motion of the mirrors during the light's round-trip 
time (quasi-static approximation). Q 

B. Free evolutions of test mass and optical field 

In this section we derive the dynamics of the free probe and the detector, i.e. that of the antisymmetric 
mode of motion of the arm-cavity mirrors when there is no light in the arm cavities, and that of the optical 
fields when the arm-cavity mirrors are held fixed. The full, coupled dynamics will be discussed in the 
following section. 

The mirror-endowed test masses are suspended from seismic isolation stacks and have free oscillation 
frequency ~ 1 Hz. However, since we are interested in frequencies above ~ 10 Hz (below these frequencies the 
seismic noise is dominant), we can approximate the antisymmetric-mode coordinate as the coordinate of a 
free particle with (reduced) mass m/4 — as is also done by KLMTV p2[[ . Hence, its free evolution is given 

by 

x< Ht)=x s + -p 8 t, (3.10) 

m 

where x s and p s are the Schrodinger operators of the canonical coordinate and momentum of the mode. 

Inserting Eq. ( 3.10 ) into Eqs. ( 2.10 ), (BTq) and using the usual commutation relations [x s ,p s ] = ih, it is 



straightforward to derive 



12 The description of a SR interferometer beyond the quasi-static approximation |33[|34| introduces nontrivial cor- 
rections to the back-action force, proportional to the power transmissivity of the input arm-cavity mirrors. Since 
the power transmissivity expected for LIGO-II is very small, we expect a small modification of our results, but an 
explicit calculation is much needed to quantify this effect. 
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As discussed in detail by BGKMTV |30), since at frequencies below < 10 Hz the data will be filtered out, 
the free evolution observable x^ ' (t) , whose Fourier component has support only at zero frequency (in a real 
interferometer it has support at the pendulum frequency ~lHz), does not contribute to the output noise. 
For this reason, henceforth, we shall disregard the free-evolution observable x^ (t) in the equations of motion 
describing the dynamics of GW interferometers. 

Concerning the free detector (the light with fixed mirrors), we can solve its dynamics by expressing the 
various quantities in terms of the quadrature operators of the input field at the SR mirror, Sj, i = 1, 2 (see 
Fig. |j) . For LIGO-II the input field will be in the vacuum state. All the quantum fluctuations affecting the 
output optical field bi are due to the vacuum fluctuations <ij entering the interferometer from the SR mirror. 



Through Eqs. (3.3), (3.8), we have already expressed Z and F in terms of the quadrature fields b^ and 
ci; thus we need now to relate the latter to Oj, i = 1,2, This can be done using Eqs. (2.11),(2.15)-(2.19) of 
Ref. |]1J|, in the case of fixed mirrors. First, for the input-output relation at the beam splitter (see Fig. [|) 
we have 

di-B ie 24/3 , d 2 =c 2 e 2t ?, (3.12) 

which is obtained from Eq. (2.11) of Ref. @, or Eq. (16) of Ref. || in the limit I Q -> and h -> 0, 
i.e. when we neglect the effects of mirror motion under radiation pressure and gravitational waves. Second, 
propagating the quadrature fields inside the SR cavity, we obtain [see Eqs. (2.16), (2.17) of Ref. H j] 

(3.13) 
(3.14) 

where </> = [iOol/c] mo d 2r: is the phase gained by the carrier frequency ujq traveling one-way in the SR cavity, 
and for simplicity we have neglected the tiny additional phase $ = VLl/ c gained by the sideband frequency 
17/27T in the SR cavity. [The length of the SR cavity is typically I ~ 10 m, hence <f> -C 1.] From the 
reflection/transmission relations at the SR mirror we derive [see Eqs. (2.18), (2.19) of Ref . Hj] 

ei = r oi + pfi , e 2 = ra 2 +ph, 

h = t fi - pax , b 2 = t / 2 - pa 2 , 

where r and p are the transmissivity and reflectivity of the SR mirror, with r 2 + p 2 = 1. 
Eqs. (3.12)-( p.l6j ) and using Eq. (3.3), we obtain for the free-evolution operators 



h 


= {d\ cos <t> — d 2 sin 4>) , 


h 


— (d\ sin (f> + d 2 cos <ji 


ei 


= (ci cos 4> + c 2 sin <j>) , 


e 2 


= (— ci sin (f> + c 2 cos 



Z { °\n)= h(tt) =^—{[(l + p 2 )cos2<f>-2pcos20\a 1 -T 2 sm2(l)a 2 } 

L J mirrors fixed iW{) 

Z ( 2 \n) = b 2 (n) = — -{t 2 sin20ai+ [(l + p 2 ) cos 20 - 2p cos 2/3] a 2 } 



[ci(n)] mi 



e 

rs fixed Mq 

[(l-pe 2 ^) cos</>Si- (1 + pe 2 ^) sin0a 2 ] 



mirrors fixed 



Mo 





(3.15) 




(3.16) 





Solving 




(3.17) 




(3.18) 




(3.19) 



13 For simplicity we ignore the effects of optical losses which were discussed in Sec. V of Ref. J14 
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where we have defined, 



M (fi) = 1 + p 2 e 4 ^ - 2p cos 20 e 2 ^ = (l + 2p cos 20 + p 2 ) ( " ° +) ( " fi - 



(ft + i 7 ) 2 



and 



Q.± = 



1 



1 + 2p cos 20 + p 2 



[±2p7 sin20-i7(l-p 2 )] 



(3.20) 



(3.21) 



Note that Z\ ' can be computed from Eqs. (3.17), ( |3.18| ) by taking the linear combination of Z\ ' and Z\ , in 



the manner of Eqs. (3^), (3__3). From Eqs. (3^) and ( 3.19 ) we obtain for the free-evolution radiation-pressure 
force: fj 



F<°>(fi) = 



2Io?iujq 



JP 



(fl 2 + 7 2 ) L 2 M 



[(l-pe 2j/3 ) COS0O! - (1 + pe 24 ' 3 ) sin0a 2 ] 



(3.22) 



Using Eqs. ( |3.17 ), ( 3. 18) ) and ( 3.22| ), and the fact that £ is frequency independent, we have explicitly checked 



that the susceptibilities of the free-evolution operators, Z\ and F^°\ satisfy the necessary and sufficient 



conditions LQM, given in Sec. II C, which define a linear quantum-measurement system with output Z. 
More specifically, using the commutation relations among the quadrature fields a^ and a 2 [Eqs. (7a), (7b) of 
Ref. pi], namely 



[ai, a 2 ,\ = — \a,2,a\,\ = 2ni S(fl — Q') , 
[oi,oj,] = = [a,!, ay] , [o2,o|,] = 



= [a 2 ,a 2 >} , 



we have derived that 



(3.23) 

(3.24) 



and we have also derived that 



R FF (n) 
R Z2F {n) = i 



2I uj 



Rz c z Q — — Rfz c 



p sin 20 



L 2 l + 2p 


cos20 + p 2 (n 

t sin0 
+ 2p cos 20 + p 

TCOS0 


-ft + )(o-n_) 

(i-p)n + »(i 
2 (n-ft+)(n- 

(i + p)n + »(i- 


•> 


. /2/ ^o 
*V hL 2 1 


-n_) 


_. _ /2I q uj q 


P)7 



TiL 2 l + 2pcos20 + p 2 (0-0 + )(fJ-Q_) ' 



(3.25) 

(3.26) 

(3.27) 
(3.28) 



i? z<F (tt) = /ZzxfCO) sinC + Rz 2 F(ty cosC- 



(3.29) 



In actuality the commutation relations ( 3.23 ), ( 3.24 ) are approximate expressions for Q <C luq- However, 
this is a good approximation since the sideband frequency 0/27T varies over the range 10 — 10 4 Hz, which is 



14 Note that if we take the limit r — *• 0, F^°'(Q) does not go to zero but ~ S(Q ± 7 tan0). Thus the main 
contribution of the fluctuating force comes from frequencies close to £7 = ±7 tan <j>, which are the optical resonances 
of the interferometer with arm-cavity mirrors fixed. 
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ten orders of magnitude smaller than ujq/2-k^ 10 14 Hz. If we had used the exact commutation relations (see 
Caves and Schumaker [||| or Eqs. (2.4), (2.5) of Ref. Q]), we would still have Rfz ( = 0, but we would 
have correction terms in the other susceptibilities. In particular, Rz ( z ( would not vanish, but would instead 
be on the order of fl/coQ. These issues are discussed in the Appendix of Ref. fl4|| . 

Before ending this section we want to discuss the resonant features of the free-evolution optical fields, which 
originally motivated the Signal Recycling (SR) j|-[|] and Resonant Sideband Extraction (RSE) schemes (j]-^. 
By definition a resonance is an infinite response to a driving force acting at a certain (complex) frequency. 
Mathematically, it corresponds to a pole of the Fourier-domain susceptibility at that (complex) frequency. 
From Eqs. ( 3.26 )— (3.29)) we deduce that Rff and Rz^f have only two poles $7±, given by Eq. (3.21), which 



are the two complex resonant frequencies of the free optical fields, Eqs. (3.17), ( |3.18 ). The corresponding 
eigenmodes are of the form e~'/ Tdcca y e _l!lo,c ', with oscillation frequency 



2p 7 sin 20 
2p cos 2<f> + p 2 



n osc± = R(n ± ) = ± - ,^L°:;7, s , (3.30) 



and decay time 



1 _ l + 2pcos2(j) + p 2 

T-docay- cj ( jj ±) - J (I - p 2 ) ' ^'^ 



This oscillation frequency and decay time give information on the frequency of perturbations to which 
the optical fields are most sensitive, and on the time these perturbations last in the interferometer before 
leaking out. Let us focus on several limiting cases: 

(i) For p — 0, i.e. the case of a conventional (LIGO-I type) of interferometer, we have ft OS c = and 
Tdecay = 1/7- Thus, there is no oscillation, while the decay time I/7 of the entire interferometer is just 
the storage time of the arm cavity. 

(ii) For p — ► 1, i.e. when the SR optical system is nearly closed, we have f2 osc = ±7tan0 and Td cca y — > +00, 
which corresponds to a pure oscillation. Noticing that for sideband fields with frequency £l/2ir, the 
phase gained in the arm cavity is 2/3 = 2 arctan fi/7 and the phase gained during a round trip in the 
SR cavity is 2<p = 2uj 1/c, we obtain that f2 osc is just the frequency at which the total round-trip phase 
in the entire cavity (arm cavity + SR cavity) is 27TO, with n an integer. 

(iii) For < p < 1 and <j) = 0, we get f2 OS c = and Tdocay = (1 + p)/[l (1 — p)\ > 1/7- This is the so-called 
tuned SR configuration [H p| , where the sideband fields remain in the intcferometer for a time longer 
than the storage time of the arm cavities [cf. (i)]. 

(iv) For < p < 1 and <j> = it/2, we get fl osc — and Td eC ay = (1 — p)/[l (1 + p)\ < 1/7- This is the 
so-called tuned RSE configuration |l]-|9|, where the sideband fields remain in the interferometer for a 
time shorter than the storage time of the arm cavities [cf. (i)]. 



It is quite straightforward to understand why Rfz c must be zero. In fact Z( is the amplitude of an outgoing 
wave; thus, the operator Z^ at an earlier time cannot be causally correlated with F at any later time, and as a 
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consequence [F (0) (ti), Zf\t 2 )] = for ti > t 2 . 



C. Coupled evolution of test mass and optical field: ponderomotive rigidity 



In Sec. [IB we have solved the equations of motion for a generic quantum-measurement device by ex- 
pressing the full- evolution operators in terms of the free-evolution operators [see Eqs. ( |2.20 )— ( 2.22|) ]. Using 
the free-evolution optical- field operators ( 3.17 ), ( 3.16 ) and ( 3.22 ) and the optical- field susceptibilities (|3.26| )- 
( 3.29 ) , along with the susceptibility of the antisymmetric mode ( 3.11 ) , qj we can now obtain the full evolution 
of the antisymmetric mode x^ 1 ) and that of the output optical field Z\ for a SR interferometer. In Ref. fl4[| , 
we evaluated the output quadrature fields by a slightly different method, introduced by KLMTV fl2j| . How- 
ever, the approach followed in this paper provides the output field in a more straightforward way, and gives 
a clearer understanding of the interferometer dynamics. Moreover, we think this method is more convenient 
when the optical configuration of the interferometer is rather complex. 

We start by investigating the interaction between the probe and the detector. The equations that couple 



the various quantities x, F and Z are [Eqs. ( 2.17)— (2.19)] : 



z { l\n) = zf\n) + R Z(F (n)x^(n) , 

F« (ft) = f<°» (ft) + R FF (Q) x (1) (ft) , 
x«(ft) = R xx (n) [G(ft) + i? (1) (ft)] . 



(3.32) 
(3.33) 
(3.34) 



In these equations, we have made explicit the dependence on the gravitational force G(ft) = — (to/4) ft 2 h(Cl) 
[see also Eq. (2.15|)] and have neglected the free evolution operator x^ ^ (see the discussion at the beginning 



of Sec. Ill B) 



Equation (3.34) is the equation of motion of the antisymmetric mode under the GW force G and the 
radiation-pressure force F, with response function R xx . Equations (3.32) and ( [3.35 ) are the equations of 
motion of the optical fields Zq and F under the modulation of the antisymmetric mode of motion of the four 
arm-cavity mirrors x, with response functions Rz c f{Q) and Rff{&), respectively. 

The optical-mechanical interaction in a conventional interferometer (p = and = 0) was analyzed by 
KLMTV in Ref. p2|. Here we summarize only the main features. Inside the arm cavity the electric field is 



[seeEq. (3.4)] 



with 



E{t) oc C cos u)ot + S\(t) cos u>ot + S% (t) sin coot . 



Sj(t) 



c 



Si(t) 

c 



cos 



2^ 



e- iQt s, 



U! t 



h.c. 



C 



3 = 1,2- 



(3.35) 



(3.36) 



where in Eq. ( 3.35 ) we have assumed that the sideband amplitudes are much smaller than the carrier 



amplitude. From Eq. (|3.35| ) we infer that the sideband fields Si and S2 modulate the amplitude and the 



As was discussed at the beginning of Sec. [II B, the free-evolution operator ar ' describing the antisymmetric 



(o) 



mode is irrelevant since it will be filtered out during the data analysis. 
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phase of the carrier field. If the arm-cavity mirrors are not moving, then it is easy to deduce that b\ ex si ex a\ 
and b 2 ex S2 ex a 2 (see Fig. |4|). Thus, given our conventions for the quadratures, we can refer to si, a\ and 
b\ as amplitude quadratures, and S2, a 2 and 62 as phase quadratures in the present case of a conventional 
interferometer. When the arm-cavity mirrors move, their motion modulates the phase of the carrier field, 
pumping part of it into the phase quadrature S 2 (t) 7 and thus into 62 [see Appendix B of Ref. ||12| , especially 
Eq. (B9a)]. As a consequence Rz 2 f 7^ but Rz x f = 0. On the other hand, the radiation-pressure force 
acting on the arm-cavity mirrors is determined by the amplitude modulation Si (t) and does not respond to 
the motion of the arm-cavity mirrors; thus Rff = 0. 

Let us now analyze a SR interferometer. As pointed out above, the antisymmetric mode of motion of the 
arm-cavity mirrors, x, only appears in the phase quadrature d 2 . [Note that now Ci and di take the place of 
ai and bi in the above analysis of conventional interferometers.] Schematically, 

avlt y i (pha 





gtlpnafle; - + U» ^ . (3.37) 



Because of the presence of the SR mirror, part of the field coming out from the beamsplitter is reflected by 
the SR mirror and fed back into the arm cavities. Due to the propagation inside the SR cavity, the outgoing 
amplitude/phase quadrature fields at the beamsplitter, di j2 , get rotated [see Eqs. ( 3.13 ), ( |3. 14 )] . Moreover, 



whereas part of the light leaks out from the SR mirror, contributing to the output field, some vacuum fields 
leak into the SR cavity from outside [see Eqs. ( 3.15| ), ( |3.16 )1. When the light reflected by the SR mirror, 



along with the vacuum fields that have leaked in, reaches the beamsplitter again, the rotation angle is 20 . 
Schematically, we can write 

cos20 — sin20\ / d± \ / vacuum fields \ (Ci\ 

sin 20 cos 20 J \d 2 ) \ from outside / \c 2 J ' 

where p and r are the amplitude reflectivity and transmissivity of the SR mirror. 

In the particular case of — or 7r/2, namely the tuned SR/RSE configurations J6|-|E)|, the rotation matrix 
in Eq. (3.38) is diagonal. Since x appears only in d 2 [see Eq. ( p.37|) ], the fact that the propagation matrix 




is diagonal guarantees that x remains only in the quadratures d 2 and c 2 - As a result, the radiation-pressure 
force, which is proportional to ci [see Eq. (|3.8| )1, is not affected by the antisymmetric mode of motion, 
and Rff = [see Eq. ( 3.26J )] as in conventional interferometers. Moreover, since the quadratures at the 



beamsplitter d\ 2 are rotated by an angle of <p when they reach the SR mirror [see Eq. ( 3.13 )1, the information 



on the motion of the arm-cavity mirrors is contained only in the output quadrature b 2 for = and b\ for 



= 7r/2. Therefore Rz x f = for = and Rz 2 f = for = w/2, as obtained directly from Eqs. (3.27) 



(3.28) 



For a generic configuration with 7^ or 7r/2, which is often referred to as the detuned case H, x appears 



in both the quadratures c±^ 2 as a consequence of the nontrivial rotation in Eq. (3.38). Thus the radiation 



pressure force and both the output quadratures respond to x, i.e. Rff 7^ and Rz i F 7^ for all £, as can 



be seen from Eqs. (3.26)-(3.2 



Before ending this section let us make some remarks. When Rff — 0, as occurs in conventional interfer- 
ometers and the tuned SR/RSE configurations, we infer from Eqs. Q3.ll), (3.33) and (|3.34) that 
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- — n 2 £ (1) (n) = G(fi) + f (0) (n) . (3.39) 

This means that the antisymmetric mode of motion of the four arm-cavity mirrors behaves as a free test mass 
subject to the GW force G(fl) and the fluctuating radiation-pressure force p(°> . It is well known that for such 
systems the Heisenberg uncertainty principle imposes a limiting noise spectral density S h = 8h/(mQ 2 L 2 ) 
for the dimensionless gravitational-wave signal h(t) = AL/L |35| . This limiting noise spectral density is called 
the standard quantum limit (SQL) for GW interferometers, and LIGO/VIRGO/GEO/TAMA interferometers 
can beat this SQL only if correlations among the optical fields are introduced |18jl9,0-nJ| . 
When R FF ^ 0, Eqs. flJUJ] ), (|3~33|) and flpg) give 



-jQ 2 x w (O) = G(0) + P (0) (O) + R FF (n) x w (fi) . (3.40) 

Thus the antisymmetric mode of motion of the four arm-cavity mirrors is not only disturbed randomly by 
the fluctuating force F^ ' , but also, and more fundamentally, is subject to a linear restoring force with a 
frequency-dependent rigidity (or "spring constant") K(Q) — —R FF (Q) ^ 0, generally called a ponderomotive 
rigidity |2fj] . This phenomenon was originally analyzed in "optical-bar" GW detectors by Braginsky, Khalili 
and colleagues, where the ponderomotive rigidity affects the internal mirror, i.e. an intra-cavity meter which 
couples the two resonators with end- mirror-endowed test masses p0[ . Hence, SR interferometers do not 
monitor the displacements of a free test mass but instead that of a test mass subject to a force field 
F res (Q) = — K(fi)x^(fi). This suggests that the SQL, derived from the monitoring of a free test mass, is 
irrelevant for detuned SR interferometers. Indeed, in Ref. JLlH we found that there exists a region of the 
parameter space p, 4> and Iq for which the quantum noise curves can beat the SQL by roughly a factor of 
two over a bandwidth Af^f. 

IV. DYNAMICS OF SIGNAL RECYCLED INTERFEROMETERS: RESONANCES AND 

INSTABILITIES 

In the previous section we have shown that in a SR interferometer the four arm-cavity mirrors are subject to 
a frequency dependent restoring force. Thus we expect the mirrors' motion be characterized by resonances 
and possible instabilities. In Refs. p3| , ^4[ , we have identified those resonances by evaluating the input- 
output relation for the quadrature fields bi (a i7 h). In this section, by using the dynamics of the whole 
system composed of the optical fields and the mirrors, we shall investigate in more detail the features of 
those resonances and instabilities. 

A. Physical origins of the two pairs of resonances 

Let us first seek a qualitative understanding of the resonances. In Fig. o we draw the amplitude and the 



phase of the ponderomotive rigidity R FF , given by Eq. ( 3.26 ), for a typical choice of LIGO-II parameters: 
(f> = 7r/2 — 0.47, p — 0.9 and Iq ~ 10 4 W. The amplitude and phase of R FF resemble those of the response 
function of a damped harmonic oscillator, except for the fact that the phase of R FF is reversed. From 
Fig. g| we infer that when the frequency / = £l/2ir is small, |-Rff| is almost constant, while the phase 
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is nearly —180°. Thus in this frequency region the spring constant is approximately a constant positive 
number ~if(Q = 0) = —Rff(Q = 0) > 0. However, K(Q, = 0) is positive only if < <j> < 7r/2, while for 
7r/2 < <j) < ir the spring constant at low frequencies is negative. As a consequence, for it/2 < (f> < tt, there 
is a non-oscillating instability, namely a pair of complex-conjugate purely imaginary resonant frequencies. 
[Note that because the SR-interferometer dynamics is invariant under the transformation — > <fi + tt |14|] , 
we can restrict ourselves to < (f> < n.] 

For larger / = 51/27T, K(fl) = —Rff{Q) has a resonant peak centered at f2 = fi osc , with width ^l/Td eC ay 
[seeEqs. ( p0| ), (§31])] • 



■2 -. 





FIG. 5. Amplitude (on the left panel) and phase (on the right panel) of Rff as a function of the sideband frequency 
f = Q./2-k for <f> = 7r/2 — 0.47, p — 0.9 and Io ~ 10 4 W. Note that the amplitude of Rff is shown in arbitrary unit. 



Hence, the dynamics of the system composed of the optical field and the arm-cavity mirrors in a SR 
interferometer is analogous to the dynamics of a massive spring, with an internal mode, attached to a test 
mass. When the test mass moves at low frequency, i.e. f2 <§; flosc, the internal configuration of the spring- 
has time to keep up with its motion and it remains uniform, providing a linear restoring force which induces 



a pair of resonances at frequencies rj mcc h = ±^4K(n <§; VL osc ) / m^ ±^J AK(£l = 0)/m. 

When the test mass moves at high frequency, the internal mode of the spring is excited, providing another 



pair of resonances to the system. Inserting the equation of motion ( 3.34 ) of x and the expression for Rff 
Eq. (3.26), into the equation of motion ( |3.33 ) of F, we obtain 

4 2Iqu>o p sin 2(f) 



(fi - n + ) (Q - 0_) F (1) (0) = driving terms - 



mtt 2 L 2 1 + 2pcos2(f> + p- 



■FW(0). 



(4.1) 



In the absence of the SR mirror, i.e. for p = 0, the term proportional to F^ on the RHS of Eq. ( j4.l| ) 
vanishes, and the optical field is characterized by the two resonant frequencies fl± given by Eq. ( 3.2l| ). By 



contrast, when the SR mirror is present, the term proportional to F^ on the RHS of Eq. (4.1) shifts the 
resonant frequencies away from the values Cl±. 

In conclusion, the dynamics of SR interferometers is characterized by two (pairs of) resonances with 
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different origin: the (pair of) resonances at iow frequency have a "mechanical" origin, coming from the 
linear restoring force due to the ponderomotive rigidity; the (pair of) resonances at higher frequency have an 
"optical" origin. Because of the motion of the arm-cavity mirrors the optical resonant frequencies get shifted 
away from the free-evolution SR resonant frequencies £l±. In this sense we can regard the SR interferometer 
as an "optical spring" [See Fig. g]. 




(1 (I (I 



V 

FIG. 6. The SR-interferometer dynamics resembles the dynamics of a massive spring with one internal oscilla- 
tion mode (and damping) attached to a test mass. The overall dynamical system is characterized by two pairs of 
resonances. 



B. Quantitative investigation of the resonances 



Equations ( |3.32 )-(3.34) describe the coupled evolution of the dynamical variables x, F and Z: 



'(«) = 



G(n) + F (0 \n) 



1 - R xx (n) R FF (n) 
p(l) M = i-R xx (n)R FF m [^) + i^)iMO)G(o) 



;(i) 



zf (n) = z™ (p.) + 



(0), 



G{Q)+F^{Q) 



(4.2) 
(4.3) 
(4.4) 



c v ; c v "" ; ' i-R xx (n)R FF (n) 

Let us first analyze these equations in the low-laser-power limit, which has long been considered in the 



literature for the SR/RSE schemes ^-§1 and has recently been tested experimentally |lfj| , pd| . For LIGO-II 
J3[ low-laser-power limit corresponds to Iq <C 10 4 W. Using Eqs. ( 3.26| )-( p.29 ), and the fact that Z\ ' does 
not depend on Iq, and F^ oc yTo [see Eqs. ( 3.17 ), ( 3.18 ) and ( 3.22J )], we deduce that R FF oc Iq and 
Rz c f oc a/^o- Therefore, for very low laser power, if we restrict ourselves only to terms up to the order of 



'To, we can reduce Eq. (4.4) to: 



zPm 



low power 



z { ,°\n) 



RzMty R xx (ty G(n) . 



(4.5) 



which says that the response of Z\ ' to the GW force G is given by the product of R xx , the response of x 
to G, times RzrFi the response of Zq to F. Hence, for low laser power the dynamics is characterized by 
four decoupled resonant frequencies: two of them, H, 2 — (degenerate), are those of the free test mass as 
embodied in R xx ; the other two, Q — Q± [see Eq. ( 3.2l| )], are those of the free-evolution optical fields as 



embodied in Rz c f- As was discussed in Sec. II B, when the imaginary part of the resonant frequency is 



2G 



negative (positive) the mode is stable (unstable). Therefore the decoupled "mechanical" resonances fi 2 = 
are marginally stable, while the decoupled "optical" resonances fi± are stable. [We remind the reader that 
S*(Ot) < 0.] 

If we increase the laser power sufficiently, the effect of the radiation pressure is no longer negligible, and 



from Eqs. (4.2)— (4.4) we derive the following condition for the resonances: 

^n(^) Rz ( f(.Q) 



1 — Rxx(Q) Rff(Q) 



(4.6) 



which simplifies to, 



fi 2 (fi-fi+)(fi-fi_) + ^p— (fi+-fi_) = o. (4.7) 

2-fsQL 

In these equations we have adopted as a reference light power i S QL = mL 2 7 4 /4w , introduced by KLMTV 
fL2| ; this is the light power at the beamsplitter needed by a conventional interferometer to reach the SQL 
at fi = 7. Because of the presence of the term proportional to Iq in Eq. (4.7), il 2 = and fi = fi± are no 
longer the resonant frequencies of the coupled SR dynamics. 



If the laser power is not very high, we expect the roots of Eq. (4.7) to differ only slightly from the decoupled 
ones. Let us then apply a perturbative analysis. Concerning the double roots fi 2 = fi 2 , = 0, working at 
leading order in the frequency shift Afio = fi — Qq = fi, we derive 

(Afi ) 2 = /o7 ' 3 ("+-"-) = _A_ (2p7 2 sin20)(l + 2pcos20 + p 2 ) 

1 °' 2/ SQL fi+fi_ /sql 4p 2 sin 2 20 + (1 - p 2 ) 2 ' l ' J 

If the SR detuning phase lies in the range < < tt/2, then (Afio) 2 is always positive. Hence, at leading 
order, the initial double zero resonant frequency fi 2 = splits into two real resonant frequencies having 
opposite signs and proportional to {Iq/Isql) 1 ^ 2 7- The imaginary parts of these resonant frequencies appear 
only at the next to leading order, and it turns out (as discussed later on in this section) that they always 
increase (becoming more positive) as Iq/Isql grows, generating instabilities. 

If the SR detuning phase lies in the range tt/2 < <f) < ir, then at leading order (Afio) 2 is negative, and we 
get two complex-conjugate purely imaginary roots. The system is therefore characterized by a non-oscillating 
instability. 



Regarding the roots fi = fi±, we can expand Eq. (4.7) with respect to Afi± — fi — fi±. A simple calculation 
gives 

/„7 3 1 



Using Eq. (3.21) we find that 



«, AO s ^ hi [4p 2 sin 2 20 - (1 - p 2 ) 2 ] (1 + 2p cos 20 + p 2 ) 2 

K(Afi ± ) = Txy r 2 . 2o , , f1 7f^r 2 , (4.10) 

2-/sql [4:p z sin 20 + (1 - p z y\ z 

5 (Att )= Io [2p7sin20(l-p 2 )](l + 2pcos20 + p 2 ) 2 
1 ±} /sql [4p 2 sin 2 20+(l-p 2 ) 2 ] 2 ■ l ' ' 

This says that, if the SR detuning phase lies in the range < < 7r/2, then 3(Afi±) always decreases 
(becoming more negative) as Iq/Isql increases. Hence, the imaginary parts of the resonant frequencies are 
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pushed away from the real f2 axis, i.e. the system remains stable. On the other hand, 3i(Af2±) may either 
increase or decrease as Io/Isql grows. If tt/2 < (j> < tt then the imaginary parts become less negative as the 
laser power increases, so the system becomes less stable. 
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FIG. 7. Shift of the resonances in a SR interferometer induced by the radiation pressure force as Iq increases from 
~0 up to Jsql- This figure is drawn for a SR mirror reflectivity p = 0.9. 

Note that, although turning up the laser power drives the optical resonant frequencies away from their 
nonzero values Q±, their changes are very small or comparable to their original values. By contrast, the 
mechanical resonant frequencies move away from zero; hence their motion is very significant. In this sense, 
as the laser power increases, the mechanical (test-mass) resonant frequencies move faster than the optical 
ones. This fact can also be understood by observing that AOo is proportional to the square root of Io, while 
Afl± is proportional to Iq itself. For the optical configurations of interest for LIGO-II, we found pi that 
when we increase the laser power from Iq = to Iq = Isql, the optical resonant frequencies stay more or 
less close to their original values while the mechanical ones, which start from zero at Iq = 0, move into the 
observation band of LIGO-II as Iq — » Isql ■ 

To get a more intuitive idea of the shift in the resonant frequencies for high laser power, we have explored 
the resonant features numerically. In Fig. |7J we plot the trajectories of the resonant frequencies when 1$ 
varies from ~ to Isql (the arrows indicate the directions of increasing power), for two choices of SR 
parameters: p = 0.9, and cf> = tt/2 =p 0.47, for which the decoupled resonant frequencies fl± coincide. The 
behaviours of the optical resonant frequencies under an increase of the power agree with the conclusion of 
the perturbative analysis deduced above. For <j) = tt/2 — 0.47, or more generally for < <j) < tt/2, the 
imaginary part of the optical resonant frequency becomes more negative when the laser power increases, and 
the resonance becomes more stable; for <f> = tt/2 — 0.47, or generically for tt/2 < <j> < tt, the imaginary part 
becomes slightly less negative when the laser power increases. The behavior of the mechanical resonance is 
particularly interesting. For <j) = tt/2 — 0.47, or generically for < <p < tt/2, and for very low laser power 
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Jo the two resonant frequencies separate along the real axis, as anticipated by the perturbative analysis. 
Moreover, as Iq increases they both gain a positive imaginary part. However, since the trajectory is tangent 
to the real axis, the growth of the imaginary parts is much smaller than the growth of the real parts. For 
<f) = 7r/2+0.47, or more generally for tt/2 < cf> < ir, the two resonant frequencies separate along the imaginary 
axis, moving in that direction as Iq increases. 

We finally note that whenever the SR detuning <fi is different from and 7r/2, the mechanical resonance is 
always unstable. We shall discuss this issue in more detail in the next section. 

C. Characterization of mechanical instabilities 

As discussed in the previous section, the coupled mechanical resonant frequencies always have a positive 
imaginary part, corresponding to an instability. The growth rate of this unstable mode is proportional to 
the positive imaginary part of the resonant frequency. The time constant, or e-folding time of the mode, is 
1/3(0). Hence, the larger the 9(0) the more unstable the system is. 

In order to quantify the consequences of the instability, we have solved numerically the condition of 



resonances, Eq. (4.7). In the left panel of Fig. K, we plot the imaginary parts of the four resonant frequencies, 



in units of 7 = Tc/AL (the bandwidth of the arm cavity, see Sec. Ill A), as a function of the detuning phase 
< < 7r of the SR cavity, fixing Iq — Isql — 10 4 W and p = 0.9. For an interferometer with arm-cavity 
length L = 4 km, and internal-mirror power reflectivity T = 0.033, which is the value anticipated by the 
LIGO-II community ||, we get 7 = 619 s _1 . Hence, the storage time of the arm cavity is 1/7 ~ 1.6 ms. 

From the left panel of Fig. H we infer that the imaginary parts of the two coupled optical resonant 
frequencies (shown with a solid line) coincide over the entire range < 4> < 7r - The imaginary parts of the 
two coupled mechanical resonant frequencies (drawn by a long-dashed line) also coincide for < 4> < 7r/2, 
but they have opposite imaginary parts for tt/2 < cf> < tt (see also Fig. for two special choices of <fi). 
From the various plots we conclude that the region characterized by the weakest instability is 4><tt/2. It is 
important to note that for these values of the detuning phase the noise curves of a SR interferometer have 
two distinct valleys that beat the SQL (see Sec. IV of Hj). Fj In Ref. |L4| the authors pointed out that 
the positions of the valleys of the noise curves coincide roughly with the real parts of the system's coupled 
mechanical and optical resonant frequencies. By taking into account Fig. |g and the dynamics of the system, 



discussed in Sec. IV A , we can make the following remark. The "spring constant" K(fl) is real only for 
O <C Cl±. For larger O's, its imaginary part contributes to that of the resonant frequency, and thus to the 
instability. Therefore, the farther the coupled mechanical resonant frequency is from the decoupled optical 
resonant frequency (0±), the less unstable it is. However, the distance between the coupled mechanical 
resonant frequency and the decoupled optical resonant frequency (0±) is directly related to the distance 
between the coupled mechanical and coupled optical resonant frequencies. Therefore, the more separate the 



17 In this paper we are only concerned with the quantum noise. Thermal noise also contributes significantly to the 
total interferometer noise; for the current baseline design it is estimated to be slightly above the SQL |15| , but design 
modifications are being explored [161 which would reduce it to about half the SQL in amplitude. 
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two coupled resonances are, i.e. the farther apart the two valleys of the noise curve are, the more stable the 
mechanical resonance is. 
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FIG. 8. The growth of instabilities for highly reflecting SR mirrors. In the left panel we plot the imaginary 
part of the resonant frequencies, obtained solving Eq. (4.7), versus the SR detuning phase <j>, for p = 0.9 and 
Io = /sql — 10 4 W. On the right panel we blow up the plot shown in the left panel for the detuning region 



V = {<j> : arctan[(47 /l! 



SQL 



,1/31 



<<j» < tt/2}, fixing p = 0.8,0.9,0.95,0.98 and I = Isql =; 10 4 W. This range of 



physical parameters corresponds to interesting LIGO-II noise curves KUJJ14J. 



In Ref. p4[ |, by analyzing the case of very highly reflecting SR mirrors (p — > 1) the authors found interesting 
noise curves for the detuning range T> = {cf> : arctan[(4Z /-7sQL) 1 ' 3 ] < <j) < tt/2} [see Sec. IV A and, in 
particular, Eq. (4.4) of Ref. 10]. In the right panel of Fig. g|, we blow up the left panel around this region 
V and plot various curves obtained by varying the SR reflectivity p = 0.8,0.9,0.95 and 0.98. We observe 
that, for this parameter set, the largest growth rate is ^0.27 ^124s _1 , corresponding to an e-folding time 
of 8 ms, which is five times larger than the arm-cavity storage time. 

Finally, we notice that the kind of instability we have found in SR interferometers has an origin similar 
to the dynamical instability induced in a detuned Fabry-Perot cavity by the radiation-pressure force acting 
on the mirrors |S||E3. 



V. CONTROL SYSTEMS FOR SIGNAL RECYCLED INTERFEROMETERS 



In this section we discuss how to suppress the instabilities present in SR interferometers by a suitable 
servo system. Since the control system must sense the mirror motion inside the observation band and act 
on (usually damp) it, there is an issue to worry about: If the dynamics is changed by the control system, 
it is not clear a priori whether the resonant dips (or at least the mechanical one which corresponds to the 
unstable resonance), which characterize the noise curves in the uncontrolled SR interferometer 13 14], will 
survive. In the following we shall show the existence of control systems that suppress the instability without 
altering the noise curves of uncontrolled interferometers, thereby relieving ourselves from the above worry. 
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FIG. 9. Scheme of the control system introduced to quench the instabilities present in a SR interferometer. The 
output Z, which contains the GW signal and the quantum noise, is sent through a linear filter with output C — Kc Z, 
and is then fed back onto the probe, i.e. the antisymmetric mode of motion of the four arm-cavity mirrors. 



A. Generic feed-back control systems: changing the dynamics without affecting the noise 



We shall identify a broad category of control systems for which, if the instability can be suppressed, the 
noise curves are not altered. We suppose that the output signal Z is sent through a linear filter Kc and then 
applied to the antisymmetric mode of the arm-cavity mirrors (see the schematic drawing in Fig. m. This 
operation corresponds to modifying the Hamiltonian (p.6|) into the form 



H = [{H v -xG) + H-d]-xF -xC . 



(5.1) 



where C is a detector observable whose free Heisenberg operator (evolving under Hrj) at time t is given, as 
required by causality, by an integration over t' < t, 



#°>(t) 



dt' K c (t-t')Z^\t'). 



(5.2) 



Physically the filter kernel Kq{t) should be a function defined for r > and should decay to zero when 
t — ► +oo. However, in order to apply Fourier analysis, we can extend its definition to r < by imposing 
Kc(t < 0) = 0, thereby obtaining 

f + OO 



C (0) (t)= / dt' K c {t-t')Z^{t'). 

J — oo 

Therefore, in the Fourier domain we have 

c i °\n) = K c {n)z {0 \n) 1 



(5.3) 



(5.4) 



where Kc(£l) is the Fourier transform of Kq{t). It is straightforward to show that the two time-domain 
properties Kq(t < 0) = and Kq(t — > +oo) — ► correspond in the Fourier domain to the requirement that 
Kc (f2) have poles only in the lower- half fi plane. 

Working in the Fourier domain and assuming that the readout scheme is homodyne detection with detec- 



tion phase C = const, we derive a set of equations of motion similar to Eqs. (3. 32)^(3.34), 
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?(0) 



z^(n) = z^ (n)+[R Z(F (Q) 



R Zc cM] 



?d) 



(fi) 



F«(n) = F^(n) + [R FF (ty + R FC( m] £ (1) (^) , 
x^(n) = R xx (n) [G(n) +fW(o) + cj 1) (f>)] , 



c™(fi) 



c£ 0) (n) + [Jfc c ,(fi) + i?c c c c («)] £ (1) («) 



(5.5) 
(5.6) 
(5.7) 
(5.8) 



Each of Eqs. (5.5), ( |5.6| ) and ( |5.8D has two response terms due to the two coupling terms between the 
probe and the detector in the total Hamiltonian (|5.l|). However, some of the responses are actually zero. In 



particular, inserting Eq. (O) into [F^(t),C^°\t')} and using the fact that [F^\t),zi°\if)] = for t > t' 



Eq. (U|)], we find Rfc ( (^) = 0. Combining Eq. Q with the fact that [zf '(i), ^ 0) (t')] = for all t, t' 



[see Eq. Q2.34D], we have R Z(C< (0) = = R Cc c c (A). Moreover, the fact that if c (t - *') = = C z co) F (o) (*, *') 
for t < t' gives the equality i?c c ir(0) = i^c(^) Rz ( f(^)- Imposing these conditions, we deduce a simplified 
set of equations of motion: 



KD, 



z ( / ! (n) 



z 



(0), 



Solving Eqs. 

x (1) (ft) = 

zf\n) = 

F«(Q) = 



. (n) + R ZcF (n)x^(n), 

FW(«) = F(°>(fi) + i? FF (ft)£ (1) (ft) , 
£«(ft) = flx»(ft) [G?(fi) + F«(ft) +Cj 1) (ft)] , 

cJ 1) (n) = if c (n)zW(n). 

[)-( [5.12 ), we obtain 

G(O) + f(°) (n) + x c (n) zj 0) (n) 

^(o) (n) + Rz c fRx X r G(n) + ~ (0) (n) 

1 — i?^ i?F_F L 

F(°)(fi) 



■Ra 



1- 


- -Rza; (^FF + Rz c F K C ) 




1 — -Rkk Rff 


1 


- Rxx (Rff + Rz ( f Kc) 




1 — Kc R xx Rz ( f 



1 — Rxx (Rff + Rz c f Kc) 



G(n) + K c z (0 \n) 



(5.9) 
(5.10) 
(5.11) 
(5.12) 



(5.13) 
(5.14) 
(5.15) 



1-KcRxxRz.F 

From the above equations ( |5.13| )-( p.l5[ ), we infer that the stability condition for the controlled system is 
determined by the positions of the roots of [1 — R X x(Rff + Rz c fKc)]- Therefore, by choosing the filter 
kernel Kq appropriately, it may be possible that all the roots have negative imaginary part, in which case 
the system will be stable. 

Before working out a specific control kernel Kc that suppresses the instability, let us notice that different 



choices of Kc give outputs ( 5.14 ) that differ only by an overall frequency-dependent normalization factor. 
This factor does not influence the interferometer's noise, since from Eq. ( 5TJ) we can see that the relative 



magnitudes of the signal (term proportional to G) and the noise (terms proportional to Z\ ' and F^ 1 ) 
depend only on the quantities inside the brackets { } and not on the factor multiplying the bracket [see 
Ref . Ej] for a detailed discussion of the noise spectral density] . Therefore if this control system can suppress 
the instability, the resulting well-behaved controlled SR interferometer will have the same noise as evaluated 
in Refs. [0,H for the uncontrolled SR interferometer. This important fact can be easily understood by 
observing that, because the whole output (the GW signal h and the noise N) is fed back onto the arm-cavity 
mirrors, h and N are suppressed in the same way by the control system, and thus their relative magnitude 
at any frequency fi is the same as if the SR interferometer had been uncontrolled. 
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B. An example of a servo system: effective damping of the test-mass 



Physically, it is quite intuitive to think of the feed-back system as a system that effectively "damps" the 
test-mass motion. When the control system is present, the equation of motion for the antisymmetric mode 



can be obtained from Eqs. (5.11), (5.9) and (5.12). It reads [as compared to Eq. (3.3 



x {1) (n)- 



ftxx 



G(ft) + F^(n) + K c z^'(a) 



(0), 



c 



1-K c R xx Rz c f 
Denoting by R%. x the response of x^ to G and F^ when the servo system is present, i.e. 

-H'XX 



R' 



c 

XX 



1 — Kc Rxx Rz c f 
we can rewrite the overall normalization factor which appears in Eqs. ( 5.13| )-(5.15) as 



R c 



(5.16) 



(5.17) 



(5.18) 



1 — Rxx (RFF + Rz c F Kc) Rxx 1 — Rxx Rff 

A sufficient condition for stability is that both R xx /R xx and 1/(1 — R xx Rpp) have poles only in the lower- 
half complex plane. [Note that when the servo system is present R xx replaces R xx in the stability condition 
of the system, see Sec. [IIB| , Eqs. ( 2.20 )-( f2.22 ) and discussions after them.] 
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FIG. 10. Effective damping due to a servo system with control kernel given by Eq. ( 5.21 ). We have fixed: A = 0.05 7, 
p — 0.9, 4> = ""/2 — 0.47 and Iq from ~ up to Isql — 10 4 W. The arrows indicate the directions of increasing light 
power Jo- The originally unstable mechanical resonance (solid line) is pushed downward in the complex Q,-plane, and 
stabilized (dashed line). The figure also shows the effect of the control system on the stable optical resonances. 
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We have found it natural to choose for R XX (Q) the susceptibility of a damped oscillator (with effective 
mass m/4), having both poles in the lower-half tt plane at fl = —iX, i.e. m 

4 1 



R xx(ty 



(5.20) 



to (n + i\) 2 ' 

with A a real parameter. This choice automatically ensures that R xx /R xx has poles only in the lower-half 
complex plane. Moreover, by choosing A appropriately we can effectively push the roots of (1 — R xx Rff) 

= 71-/2 



in Eq. (5.18) to the lower-half O plane, as shown in Fig. |10| for p = 0.9, 
from — up to /sql- 



0.47, A = 0.05 7 and I 



However, we also need to check that Kc(£l) has poles only in the lower-half Q plane. Using Eqs. (5.17) 



(5.20) we obtain the following explicit expression for the kernel: 



K c (fl) 




(i + 2 P cos2(j) + P 2 )(n ~ n_)(o - n+) 

(ft + ij) cos(0 + C) + P(M - il) cos(0 - C) 



(5.21) 



For C, — or C, = w/2, i.e. when either of the two quadratures b\ or bi is measured, the control kernel ( 5.21 ) 
indeed has poles only in the lower-half complex plane. More generally, we have shown that if < <f) < t/2, 



the control kernel (5.21) has poles in the lower-half complex plane for all n/2 < £ < n, regardless of the 
value of p, but it may become unphysical in the region < £ < tt/2. However, for the unphysical values of 
£ there are various feasible ways out. For example, we could change R xx by replacing to in Eq. ( 5. 20] ) with 
a slightly smaller quantity mc- In this case 



R.i 



1 
W 



TO 

T 



n i 



mc \ 

TO / 



iX 



n i 



+ iX 



(5.22) 



By choosing mc appropriately, we can use the first factor in Eq. ( 5.22J ), which has a root in the upper-half 



complex plane, to cancel the bad pole coming from Rz ( f in Eq. ( 5.21 ), so that K~c will have poles only in the 
lower-half complex plane. Finally, we must adjust A so that the effective damping suppress the instability. 

Of course, the servo electronics employed to implement the control system will inevitably introduce some 
noise into the interferometer. In our investigation we have not modelled this noise. However, LIGO exper- 
imentalists have seen no fundamental noise limit in implementing control kernels of the kind we discussed, 
and deem it technically possible to suppress any contribution coming from the electronics to within 10% of 
the total predicted quantum noise |37| , |38| . This issue deserves a more careful study and it will be tackled 
elsewhere p9| . 

In this paper we have restricted ourselves to the readout scheme of frequency independent homodyne 
detection, in which only one (frequency independent) quadrature b^ is measured. The issue of control-system 



In the time domain this choice of R^. x (fi) corresponds to the equation of motion 



m .. 



m\ 



mX 



X x + forces . 



(5.19) 



:-S4 



design when other readout schemes are present, e.g., the so-called radio- frequency modulation-demodulation 



design, is currently under investigation 39 



Finally, for simplicity we have limited our discussion to lossless SR interferometers. When optical losses 
are taken into account, we have found that the instability problem is still present |1J] and we have checked 
that those instabilities can be cured by the same type of control system as was discussed above for lossless 
SR interferometers. 

VI. CONCLUSIONS 

Using the formalism of linear quantum-measurement theory, extended by Braginsky and Khalili pM to 
GW detectors, we have described the optical-mechanical dynamics of SR interferometers such as LIGO-II H. 
This analysis has allowed us to work out various significant features of such interferometers, which previous 
investigations pHq] could not reveal. 

We have found that when the (carrier) laser frequency is detuned in the SR cavity, the arm-cavity mirrors 
are not only perturbed by a random fluctuating force but are also subject to a linear restoring force with a 
specific frequency-dependent rigidity. This phenomenon is not unique to SR interferometers; it is a generic 
feature of detuned cavities |36|,|2(],|33],[34| and was originally used by Braginsky, Khalili and colleagues in 
designing the "optical bar" GW detectors |2(| ■ 

Our analysis has revealed that, for SR interferometers, the dynamics of the whole optical-mechanical 
system, composed of the arm-cavity mirrors and the optical field, resembles that of a free test mass (mirror 
motion) connected to a massive spring (optical fields). When the test mass and the spring are not connected 
(e.g., for very low laser power) they have their own eigenmodes, namely the uniform translation mode for 
the free test mass (free antisymmetric mode), and the longitudinal- wave mode for the spring (decoupled 
SR optical resonance). However, as soon as the free test mass is connected to the massive spring (e.g, for 
LIGO-II laser power), the two free modes get shifted in frequency, so the entire coupled system can resonate 
at two pairs of finite frequencies (coupled mechanical and optical resonances) . From this point of view a SR 
interferometer behaves like an "optical spring" detector. For LIGO-II parameters, both resonant frequencies 
can lie in the observation band 10 Hz < / < 10 kHz and they are responsible for the beating of the SQL in 
SR interferometers ||Q. 

The formalism used in the present paper has allowed us to analyze in more detail the features of the insta- 
bilities in SR interferometers, pointed out in Refs. 0,0|. Most importantly, we have shown the possibility 
of using a feed-back control system to cure such instabilities without compromising the performance of the 
interferometer. However, before any practical implementation, a much more careful and precise study should 
be carried out, including various readout schemes p9| . 

Finally, the general discussion based on the Braginsky-Khalili force-susceptibility formalism, given in the 
first part of this paper (Sec. O), and the application to a specific type of GW interferometer, the LIGO-II 
SR interferometer, given in the second part of it (Sees. |IH-M), may provide, along with Refs. |0,n, a 
framework for future investigations of quantum noise in advanced, more complex, optical configurations. 
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APPENDIX A: BASIC PROPERTIES OF LINEAR SYSTEMS 

In this Appendix, to clarify the formalism used in Sec. 0, we summarize some well-known basic properties 
of linear systems linearly coupled to each other or to external classical forces. Much of this material can be 
found in Sakurai [J40|, and for its application to quantum-measurement processes in Braginsky and Khalili 
|2l| and Caves et al. ||. 

Definition 1 (Linear systems) Any system whose Hamiltonian is at most quadratic in its canonical co- 
ordinates and momenta is a linear system. 



Definition 2 (Linear observables) Any linear combination (either time dependent or time independent) 
of the canonical coordinates and momenta of a linear system, plus a possible complex number (C-number), 
is a linear observable of the system. 

Denoting all the canonical coordinates and momenta by Ci with i = 1,2, • • •, the Hamiltonian of a linear 
system can be written as 



H(t) = J2 L 2( t )^C J +J2L\(t)C l + L (t) . 



(Al) 



where L^ W i s symmetric in i and j. The equations of motion of the canonical observables in the Hcisenberg 
picture read [we use the fact that Cjh does not depend explicitly on time]: 



ih—C jH (t)= \C jH (t),Hff(t) 



i/t(-oo,t) [c jS ,H s (tj\ U(-oo,t), 



W(-oo,t) 



^ 2 4 m (i)C^C mg +$>!(*) Q 



U(-oo,t), 



= J2 2 L 2 m (t) Cji C mH (t) + Y, L i(t) d j 



(A2) 



/ . m 



Here the subscripts S and H stand for Schrodinger and Heisenberg pictures respectively, Cji = [Cjs,Cis] is 
the commutator between the canonical operators, which is a C-number, and U(— oo, t) is the time-evolution 
operator which satisfies the Schrodinger equation 



ih—U(-oo, t) = H s U(-oo, t) 



with initial condition U{— oo, — c») = 1. The solution to Eq. ( |A2[) is of the form 

CjH (t) = Yl a Jk (*) Cuh (-oo) + I3 3 (t) = Y Oijk (t) C kS + 0j (t) , 



(A3) 



(A4) 



where ctjk{t) and (3j(t) are time dependent C-numbers. 

For any linear observable A it follows from linearity that Ajj(t) = J2i a j(t)CjH(t) + b(t), which, along 



with Eq. (A.4), leads to 
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Ah (t) = J2 a J (*) Cm (*) + 5 W = Z! a J W Q J fc (*) ^s + J2 a i (*)#»■ (*) + ft (*) • 

3 J,* J 



(A5) 



This provides the following theorem: 



Theorem 1 At any time the operator of a linear observable in the Heisenberg picture can always be written 
as a linear combination of operators of the (time-independent) canonical variables in the Schrodinger picture 
plus a possible C-number. 

Applying the above theorem to any two linear observables A and B, recalling that Cjk = [Cjs,Cks] is a 
C-number and the commutator between a C-number and any operator is zero, we find 

A H (t),B H (t')] = J2 if H) lk(t') C 3k , (A6) 

which is a C-number. Therefore, the following theorem holds: 

Theorem 2 In the Heisenberg picture, the commutator of the operators of any two linear observables at two 
times is a C-number. 

We are interested in the evolution of a linear system subject to a classical external linear force or linearly 
coupled to another independent linear system. A force-susceptibility kind of formulation can be introduced in 
these cases (as is done by Braginsky and Khalili, see Sec. 6.4 of Ref. pifl ). We shall describe the system using a 
perturbative approach. Thus we write the total Hamiltonian in the Schrodinger picture as H$ = Hos + Vs(t), 
where Vs(t) is treated as a perturbation with respect to the zeroth order Hamiltonian Hqs- It is generally 
convenient to introduce the so-called Interaction picture (see, e.g., Sections 5.5 and 5.6 of Ref. fhj]), in which 
the evolution operator Uj is defined by the relation U(— 00, t) = Uo(—oo,t) Uj(— 00, i), where Uq(— 00, t) is 
the evolution operator associated with H s and U is defined by Eq. <\K3j ). Then, Ui(— 00, t) satisfies the 
equations 



ih— J7i(-oo,«) = Vi(t) Ui(-oo, t) , 
at 



Uj(— 00, —00) = 1 



(A7) 



with Vr(t) = Uq (— 00, t) Vs(t) Uo(—oo, t). The solution of Eq. (A7) can be written as a perturbative expan- 
sion, 



Uj{-oo,t) = 1 + - 



= ±-(- 

^-^ tj.I \ if). 



n=0 



dhVjih) 



T ■ 



dti 



dt 2 V I {t x )V I {t 2 



dhV^h 



(A8) 



where T denotes the time-ordered product H] . The Heisenberg operator associated with any observable A, 
evolving under the full Hamiltonian H, is linked to the corresponding Heisenberg operator evolving under 
the Hamiltonian Hq by the relation Aji(t) = Uj(— 00, t)Ajj'(t)Ui(— 00, t), where the superscript (0) on the 
observable A denotes that the evolution is due to H . Inserting Eq. (|A8[) into the above equation, we get 



(0), 



A H (t) = A^(t) + 



dti 



dti 



V^hlA^it) 

■tn-l 



dti 



dti 



(0), 



Vj(t 2 ), Vzit^A™® 



dti-- 



dt, 



Vi(t n ) 



Vjih), Vi(t!),A%>(t) 



(0), 



(A9) 
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For a linear system subject to an external classical linear force G(t), the interaction term is Vi(t) — —x H G(t). 



Plugging this expression into Eq. (A9) and using Theorem 2, it is straightforward to deduce that the second 
and all higher order terms in Eq. flA9| ) vanish and the first order perturbation gives the exact solution. 
Hence, we obtain the following theorem: 

Theorem 3 Consider a linear system subject to a classical generalized force G(t), whose Hamiltonian is 
given by H = H — x G(t), where x is a linear observable. Then, for any linear observable A, the Heisenberg 
operator An(t) can be written as the sum of its free- evolution part, A H (t), plus a term which is due to the 
presence of the external force, i. e. 

A H (t) = Af{t)+ l -J dt'C Ax {t,t')G{t'), (A10) 

where C Ax {t,t') is a C-number, called the (time-domain) susceptibility, given explicitly by 

C Ax (t,t') = [A%\t),x { °\t')}. (All) 

Let us now suppose that we have two independent linear systems V (e.g., the probe) and V (e.g., the 
detector), which by definition are described by two different Hilbert spaces Tip and Jip,. We introduce the 
Hilbert space Ti = Tip £g> Hd and define for any operator x of the system V the corresponding operator acting 
on Ti as x (g> 1 , while for any operator F of the system T> we introduce the operator 1 £g) F which acts on Ti. 
Henceforth, we shall limit ourselves to interaction terms V, in the total Hamiltonian H = Hp + H-p + V, 



of the form: V = — x ® F, with x and F acting on V and V, respectively. Using Eq. (A9) with Viit) — 
~x H (t)Ffl (t), noticing that (i) the zeroth order Heisenberg operators of two observables living in different 
Hilbert spaces commute and (ii) the zeroth order Heisenberg operators of two linear observables living in 
the same Hilbert space have a C-number commutator, we derive the following theorem: 

Theorem 4 Consider two independent linear systems V and T>, and two linear observables, x of V and F 
of T> . Suppose that the two systems are coupled by a term —x® F, i.e. the Hamiltonian of the composite 
system V + T> reads H = H-p + H-p — X®F. Then, for any linear observable A of the system V and B of 
the system T>, their full Heisenberg evolutions are given by: 

A H {t) = A%\t)+ l - J dt'C Ax (t,t')F H (t>), B H {t) = B ( °\t)+ l - J dt'C BF (t,t')x H (t'), (A12) 

where A H and B H stand for the free Heisenberg evolutions, and the susceptibilities are defined by 

C Ax {t,t') = [A%\t),x%\t')] , C BF (t,t') = [B%\t),F<®{t)] . (A13) 

In the case where the zeroth order Hamiltonian is time independent, it is easy and convenient to express the 
above formalism in the Fourier domain. We first notice that for a time independent Hq, Uo(t, t+r) = er lHaT l h 
and for any two linear observables A\ and A^ we have C AlA2 (t + r,t' + t) = C AlA2 (t,t'), i.e. C AlA2 (t,t') 
depends only on t — t' . Defining the Fourier transform of any observable A(t) as 

/•+oo 

A{Q)= / dte mt A(t), (A14) 



:W 



Eq. QA1C| ) becomes A H (tt) = Afi'(£l) + R Ax (ty G(il) while Eq. QA12| ) can be recast in the form 



\i°)fn\ , u. fn\G..fn\ a../m_S(o) 



A H (n) = A^(a) + RA X (n)F H (n), B H (n) = B [ H "{n) + R BF (n)s H {n), (ais) 



where R AB (£l) is the susceptibility in the Fourier-domain, given by 

. i^+oo ■ ,.+00 

Rab({1) = ~ / dTe mT Q(T)C AB (0,-T) = - / dre mr C A 

" J-oo "> JO 



(0,-r), (A16) 



with 0(t) the step function. For future reference, let us point out two properties which R AB (Q) satisfies 
and that we use repeatedly in Sec. ||: 

R* AB (Q.) = R AB (-Cl) , [^(niJ.Bg^fia)] =-2mh6(Q 1 +n 2 )[R AB (Q 1 )-R BA (n 2 )]. (A17) 



To deduce the first identity in Eq. (A.17), we consider the complex (Hcrmitian) conjugate of Eq. (A16) 
and use the Hermiticy of AS H '(t) and Bg(t). For the second identity in Eq. (A17), we take the double 
Fourier transform of [Ajj (t\), Bjjfa)] with respect to t\ and t%, and then using Eq. ( A16 ) we find that the 



region corresponding to t\ > t 2 in the double integral yields the R AB term of Eq. (A17), while the region 
corresponding to t\ < t 2 gives the R BA term. 
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